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Objectives 


The  objectives  of  this  work  were  to  measure  reactive  chlorine  (Cl,  CIO,  and  C1202)  and 
carbon  dioxide  (CO2)  and  to  examine  the  mechanisms  of  ozone  loss  in  the  plumes  of 
various  rockets  as  part  of  an  ongoing  investigation  of  the  impacts  of  rocket  emissions  on 
stratospheric  ozone.  Following  a  budget  revision  and  change  of  priorities  of  the  mission 
scientist  during  the  RISO  campaign,  the  project  was  rescoped  to  emphasize  ultra-fast 
response  measurements  of  Cl  and  CIO  with  ancillary  measurements  of  C02  in  rocket 
plumes,  with  some  ambient  sampling  to  support  the  NAS  A/NO  AA  ACCENT  campaigns 
of  1999  and  2000. 


Accomplishments/New  Findings 

At  the  close  of  this  3-year  project,  we  successfully  met  all  of  the  objectives  outlined  in 
the  original  proposal,  having  demonstrated  that  ozone  losses  in  rocket  plumes  proceed  at 
rates  that  can  be  explained  by  ’standard’  ozone  loss  chemistry  involving  theCIO  radical  - 
chemistry  that  is  analogous  to  that  occurring  in  the  winter  polar  vortices.  This  work  has 
also  helped  to  identify  an  error  in  a  key  kinetic  parameter  (the  rate  constant  for  the  C10+ 
reaction  CIO)  as  first  reported  by  Ross  et  al.  (2000).  In  combination  with  fast  response 
measurements  of  Cl  and  CIO  in  the  plume  of  STS- 106,  we  have  shown  that  this  rate 
parameter,  along  with  consideration  of  non-linear  effects  can  explain  ozone  loss  rates 
within  the  uncertainties  of  the  measurements  (unpublished  results  described  in  Year  4 
Annual  Progress  Report). 

An  unexpected  result  from  this  effort  was  the  discovery  of  volatile  particles  likely 
composed  of  HNO3  and  H20  in  the  fresh  plumes  of  several  small  rockets  (Atlas  HAS  and 
Athena  II).  These  particles  persisted  long  after  the  plume  had  diluted,  indicating  that  our 
understanding  of  evaporative  properties  of  binary  (or  ternary?)  particles  is  lacking.  These 
results  have  important  implications  for  evolution  of  the  plume,  including  visibility  of 
plumes  and  chemistry  both  in  the  plumes  and  as  they  mix  with  background  air.  To  date, 
these  potential  effects  are  not  included  in  models  of  plume  impacts. 

These  results  are  now  presented  in  5  publications  that  have  appeared  in  the  peer-reviewed 
literature,  four  of  which  have  appeared  in  press  and  are  attached  to  this  file  for  brevity. 
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[i]  Carbon  dioxide  (C02)  and  large  particles  (0.37-4.0  pm)  were  measured  in  the 
stratospheric  plume  wakes  of  three  rockets,  an  Atlas  HAS,  a  Delta  II,  and  an  Athena  IL 
The  correlations  between  C02  mass  and  particle  number  densities  in  each  plume  are 
consistent  with  the  unique  combination  of  solid  and  liquid  engine  emissions  of  each 
rocket  Measured  size  distributions  indicate  a  LI  pm  mode  with  density  of  2  g  cm“3, 
consistent  with  spherical  alumina  particles  emitted  by  solid  rocket  motors.  Disagreement 
between  the  measured  size  distributions  and  the  mean  sizes  inferred  from  the  known 
alumina  and  C02  emission  indices  and  an  observed  increase  in  the  particle  number 
emission  index  with  altitude  are  evidence  for  large  particle  oversampling  effects  and  the 
presence  of  condensed  volatile  compounds  within  the  particle  population.  Direct  evidence 
for  the  latter  is  a  persistent  ~0.5 — 1  part  per  million  (ppm)  shortfall  of  water  vapor  relative 
to  C02  measured  in  the  plume  of  the  Athena  II  rocket  based  on  the  expected  H20/C02 
emission  ratio.  Although  pure  ice  particles  would  not  persist  at  the  conditions  of  the 
measurements,  a  more  stable  coating  of  HN03  (as  either  nitric  acid  trihydrate  or  as  a  liquid 
layer)  could  have  reduced  the  sublimation  rate  of  the  underlying  ice,  thereby  increasing  the 
lifetimes  of  volatile  particles  within  the  plume.  If  confirmed,  such  a  process  would  have 
important  implications  for  the  radiative  and  chemical  properties  of  rocket  plumes, 
including  global  ozone  depletion  associated  with  rocket  launch  activities.  index  terms : 

0305  Atmospheric  Composition  and  Structure:  Aerosols  and  particles  (0345,  4801);  0394  Atmospheric 
Composition  and  Structure:  Instruments  and  techniques;  0340  Atmospheric  Composition  and  Structure: 

Middle  atmosphere-composition  and  chemistry;  KEYWORDS:  carbon  dioxide  (C02),  H20,  ozone  depletion, 
rocket  plumes,  emissions 

Citation:  Gates,  A.  M.,  et  ai.  In  situ  measurements  of  carbon  dioxide,  0.37  -4.0  pm  particles,  and  water  vapor  in  the  stratospheric 
plumes  of  small  rockets,  J  Geophys .  Res.,  I07( D22),  4649,  doi :  10.1 029/2GO2JDG02 121,  2002. 


1.  Introduction 

[2]  Rockets  used  in  the  space  transportation  Industiy  emit 
a  variety  of  exhaust  products  directly  into  the  middle  and 
upper  atmosphere.  The  potential  effects  of  these  emissions 
on  climate  and  stratospheric  ozone  abundances  were  first 
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discussed  during  the  1975  Climate  Impact  Assessment 
Program  [see  Hoshizaki,  1975],  and  there  has  been  subse¬ 
quent  progress  in  understanding  the  composition  of  rocket 
engine  emissions  and  their  impacts  on  the  stratosphere. 
Various  models  of  the  global  atmosphere  are  in  broad 
agreement  that  present-day  global-averaged  ozone  loss  from 
rockets  is  a  few  hundredths  of  a  percent,  a  small  fraction  of 
that  believed  to  be  due  to  all  industrial  activities  [World 
Meteorological  Organization  ( WMO ),  1992].  However,  the 
representation  of  rocket  emissions  by  the  models  has 
focused  primarily  on  the  solid  rocket  motors  (SRMs)  while 
excluding  liquid  propellants  [Prather  et  al ,  1990,  Jackman 
et  al,  1998;  Danilin  et  al ,  2001],  Furthermore,  the  models 
may  not  include  all  of  the  important  interactions  between 
the  rocket  emissions  and  natural  stratospheric  processes. 
Therefore,  until  models  are  validated  by  observations 
directly  within  the  stratospheric  plumes  of  rockets,  their 
predictions  of  ozone  losses  should  be  considered  as  lower 
limits  to  those  that  could  be  expected. 

[3]  Based  on  model  predictions  of  the  stratospheric 
impacts  of  SRM  emissions,  in  particular  their  small  poten¬ 
tial  to  deplete  ozone  relative  to  other  human  activities 
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[ WMO ,  1992],  rocket  engine  emissions  have  remained 
unregulated  [United  Nations  Environment  Programme , 
2001],  This  status  relies  upon  the  validity  of  a  number  of 
assumptions,  including  (1)  that  global  models  accurately 
predict  ozone  losses  due  to  SRM  emissions,  (2)  that  ozone 
losses  from  liquid  propellants  do  not  substantially  exceed 
those  from  SRMs,  and  (3)  that  future  rocket  emissions 
remain  comparable  to  those  at  present.  The  commercial 
space  transportation  industry  has  grown  much  faster  than 
the  economy  as  a  whole  in  recent  years  and  is  projected  to 
continue  to  do  so.  Meanwhile,  new*  liquid  propellant  rockets 
are  being  developed.  Therefore,  the  Impact  of  rocket  emis¬ 
sions  on  stratospheric  ozone  must  be  more  carefully  exam¬ 
ined  in  order  to  continue  to  meet  the  evolving  information 
needs  of  policy  makers  [Ko  et  al ,  1994]. 

[4]  To  date,  efforts  to  model  the  global  impacts  of  rockets 
have  focused  on  the  gas  (HC1  and  CU)  and  particulate 
(alumina,  or  A1203)  emissions  from  SRMs  that  reduce 
stratospheric  ozone  by  enhancing  abundances  of  inorganic 
chlorine  and  increasing  the  reactive  surface  area  of  particles 
[Prather  et  at ,  1990,  Jackman  et  al,  1998;  Danilin  et  al , 
2001],  Chlorine  chemistry  has  been  studied  extensively  in 
recent  years  in  an  effort  to  assess  the  impact  of  chlorofluor- 
ocarbons  (CFCs),  such  that  the  regional  and  global  influ¬ 
ences  of  gaseous  chlorine  compounds  emitted  by  SRMs  are 
reasonably  well  understood.  Consequently,  there  is  general 
agreement  as  to  the  extent  of  global  stratospheric  ozone 
losses  associated  with  the  increases  of  inorganic  chlorine 
abundances  due  to  SRMs.  In  contrast,  the  global  impact  of 
particles  emitted  by  SRMs  is  uncertain  and  may  even  be 
more  significant  than  that  due  to  the  gaseous  chlorine 
emissions  [Danilin  et  al,  2001], 

[5]  A  variety  of  heterogeneous  reactions  have  been  pro¬ 
posed  that  promote  conversion  of  reservoir  forms  of  inor¬ 
ganic  chlorine  into  ozone  destroying  forms  on  alumina 
particles,  for  example. 


CINO3  +  HC1  —  CI2  +  HN03t 

(I) 

CINO3  +  H20  —  HOC1  +  HN03, 

(2) 

HCl  +  HOCl->Cl2-FH20, 

(3) 

though  the  list  is  not  exhaustive,  and  only  reaction  (1)  has 
been  studied  in  any  detail  in  the  laboratory  [Molina  et  al , 
1997].  Recent  modeling  studies  [Jackman  et  al,  1998; 
Danilin  et  al,  2001]  have  demonstrated  that  ozone  loss  is 
critically  sensitive  to  the  fraction  of  alumina  mass  contained 
in  particles  small  enough  (<0.4  pm)  to  have  atmospheric 
residence  times  exceeding  a  few  months.  In  these  models, 
when  the  mass  fraction  of  submicron  particles  accounts  for 
more  than  a  few  percent  of  the  total  alumina  mass,  the 
predicted  ozone  losses  associated  with  heterogeneous 
reactions  on  alumina  exceed  those  due  solely  to  the  increase 
in  chlorine  abundances.  Although  size  distributions  of 
alumina  particles  emitted  by  rockets  have  been  studied  by 
various  researchers  [Brady  and  Martin,  1995;  Cofer  et  al, 
1987;  Dao  et  al,  1997;  Belting ,  1997;  Ross  et  al,  1999], 
there  is  no  consensus  on  the  relative  mass  fraction  of  the 
submicron  portion  of  the  size  distribution.  Reported  values 
for  these  important  particles  have  ranged  from  0.01%  to 
95%  of  the  mass  of  emitted  alumina.  Consequently,  lacking 


a  clear  understanding  of  the  alumina  size  distribution,  the 
estimates  of  the  impact  of  SRMs  are  highly  uncertain, 

[5]  Additional  uncertainty  concerns  the  chemical  nature 
of  the  alumina  surface  at  stratospheric  conditions.  The 
reactivity  of  alumina  particles  depends  on  the  nature  of 
their  surfaces,  Including  the  nature  and  availability  of 
condensed  materials.  Limited  laboratory  studies  of  alumina 
reactivity  in  the  presence  of  HC1  and  H20  [Molina  et  al, 
1997;  Nelson  et  al,  2001]  suggest  that  alumina  will  adsorb 
volatile  compounds  which  then  determine  the  kinetics  of  the 
various  heterogeneous  reactions.  Such  volatiles  could  be 
derived  from  the  stratospheric  background  or  from  com¬ 
pounds  emitted  by  the  rocket  itself.  Because  condensed 
phases  are  known  to  play  critical  roles  in  stratospheric 
chemistry,  it  is  Important  to  understand  the  nature  of  the 
interactions  of  volatile  compounds  with  alumina  particles, 
[7]  In  this  paper  we  report  on  the  first  simultaneous 
measurements  of  particles  and  CO2  In  the  plume  wakes  of 
three  rockets  that  employ  different  mixtures  of  solid  and 
liquid  propellants.  Previous  measurements  of  SRM  alumina 
particles  were  limited  by  the  lack  of  a  conservative  tracer  of 
the  rate  of  propellant  consumption,  which  precluded  vali¬ 
dation  of  the  observed  particle  size  distributions.  The 
absolute  emission  of  C02  is  well  understood  In  solid  and 
liquid  propellant  rocket  engines,  such  that  this  species 
serves  as  an  excellent  tracer  to  understand  the  emissions 
of  less  well-understood  rocket  combustion  products  and 
their  impacts  [Schulte  and  Schlager,  1996].  We  use  calcu¬ 
lations  from  an  after-burning  rocket  plume  flow  field  model 
to  establish  expected  ratios  of  the  deposition  rates  of  various 
plume  species  to  that  of  C02  [e.g,,  see  Zittel ,  1994;  Denison 
et  al,  1994].  Previous  measurements  also  were  not  able  to 
distinguish  between  nonvolatile  (alumina)  and  volatile  (e.g., 
ice  compounds)  components  of  the  plume  particle  popula¬ 
tion.  Here  we  examine  correlations  between  C02  and 
particles  and  between  C02  and  H20,  as  well  as  particle 
size  distributions,  to  infer  some  of  the  characteristics  of 
alumina  in  SRM  plumes.  We  discuss  some  important 
implications  that  these  characteristics  could  have  for  assess¬ 
ments  of  the  Impact  of  rocket  emissions  on  stratospheric 
ozone. 


2,  Observations 

[s]  The  Atmospheric  Chemistry  of  Combustion  Emis¬ 
sions  Near  the  Tropopause  (ACCENT)  campaigns  took 
place  In  April  and  September  1999,  during  which  the 
NASA  WB-57F  high-altitude  aircraft  carried  Instruments 
to  measure  a  wide  variety  of  gas  and  aerosol  species  in 
rocket  plume  wakes.  In  this  paper  we  limit  our  attention  to 
measurements  of  C02,  particles  with  diameters  between 
0,37  and  4,0  pm,  water  vapor,  and  ambient  pressures  and 
temperatures, 

[9]  Data  were  obtained  In  the  plume  wakes  of  Atlas  HAS, 
Delta  II,  and  Athena  II  launch  vehicles  that  use  different 
combinations  of  solid  and  liquid  (kerosene  and  liquid  oxy¬ 
gen)  propellant  rocket  engines,  thereby  producing  different 
relative  emissions  of  C02,  H20,  and  alumina,  as  summar¬ 
ized  in  Table  1.  Atlas  HAS  plume  data  were  obtained 
between  22,97  and  23,92  UT  (solar  zenith  angles,  or 
SZA,  >87°)  on  April  12,  1999,  following  launch  at  22.83 
UT  from  Cape  Canaveral,  Florida  (28°33/N,  80°  18'W). 
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Table  1.  Emission  Mass  Ratios  Relative  to  C02 


Rocket 

Altitude, 

km 

ai2o3, 
gfg  co2 

h2o, 

g/g  CO, 

HC1, 
g/g  co. 

Athena  II 

16-19 

0,94 

0.78 

0.36 

Atlas  HAS 

19 

0.17 

0.47 

0.06 

Delta  II 

19 

0.37 

0,57 

0,12 

11-13 

0,55 

0.65 

0.24 

Delta  II  plume  data  were  obtained  between  18.61  and  19,41 
UT  (SZA  —34°)  following  launch  at  1 8,53  on  April  1 5}  1999, 
from  Vandenberg  Air  Force  Base,  CA  (34°  4814, 120°37/W). 
Athena  II  plume  data  were  obtained  between  1 8,43  and  1 8,98 
UT  (SZA  —42°)  following  launch  from  Vandenberg  AFB  at 
18,35  UT  on  September  24, 1999,  Similar  to  previous  studies 
l Ross  etal ,  1997, 2000],  die  WB-57F  (airspeed  -180  m  s"1) 
flew  though  each  plume  wake  a  number  of  times  with 
encounters  ranging  from  10-70  s  in  duration.  In  general, 
each  encounter  was  flown  through  different  portions  of  the 
plume  wake,  and  we  assume  that  gas  and  particle  enhance¬ 
ments  are  due  solely  to  rocket  engines,  rather  than  to  engine 
exhaust  from  the  WB-57P. 

[io]  The  Atlas  HAS  plume  was  encountered  fourteen 
times  near  19  km  altitude,  and  the  WB-57F  aircrew  reported 
excellent  plume  visibility  for  about  60  min  until  sunset; 
however,  a  complete  set  of  measurements  was  obtained  only 
during  two  encounters.  The  Delta  II  plume  was  encountered 
12  times,  the  first  occurring  at  an  altitude  of  19  km.  After 
the  first  encounter,  the  plume  was  no  longer  visible  to  the 
WB-57F  aircrew  at  this  altitude  and  the  aircraft  descended 


to  the  11-12  km  region  where  the  plume  maintained 
visibility.  At  these  lower  altitudes  the  Delta  II  bums  six 
SRMs  along  with  the  liquid-fueled  core  stage,  a  configu¬ 
ration  that  differs  from  that  at  19  km,  where  only  three 
SRMs  bum  with  the  core  stage.  The  Athena  II  plume  was 
encountered  six  times,  the  first  five  near  19  km.  After  about 
20  min,  the  aircrew  reported  difficulty  locating  the  plume, 
and  the  aircraft  descended  to  17  km  for  one  last  encounter 
before  landing.  Following  a  discussion  of  the  nature  of  the 
observed  partiele-C02  and  H20-C02  relationships,  we  offer 
an  explanation  for  these  high-altitude  visibility  issues. 
Important  quantities  for  all  of  the  plume  encounters  are 
summarized  in  Table  2. 

3.  Measurements 

[u]  C02  abundances  were  measured  at  5  Hz  by  infrared 
absorption  using  a  nondlspersive  infrared  detector  (Li-Cor, 
Model  6251,  Lincoln,  Nebraska)  repackaged  for  the  low- 
temperature  (—195-220  K)  and  low-pressure  (60-200  hPa) 
WB-57F  aircraft  environment.  The  instrument  detects  dif¬ 
ferences  in  concentrations  between  two  cells  held  at  the 
same  temperature  (— 30-40°C)  and  pressure  (~400  hPa). 
Under  controlled  laboratory  conditions,  this  instrument 
detects  changes  in  C02  mixing  ratios  of  less  than  0.1  ppm 
in  less  than  I  s.  On  an  airborne  platform,  measuring  with 
this  sensitivity  requires  extremely  precise  pressure  and 
temperature  control.  During  flight  a  constant  mixing  ratio 
is  maintained  in  one  of  the  cells  (the  “reference”)  by 
flowing  a  C02  standard  of  known  concentration  near  365 


Table  2*  Details  of  Plume  Encounters 


A 

Start, 
UT  sec 

B 

End, 
UT  sec 

C 

T, 

K 

D 

P» 

hPa 

E 

mi 

1018  cm”3 

F 

Theta, 

K 

G 

Mas  C02, 
_ PPm 

H 

Int  CO,, 

10“8  g  s  cm”3 

I 

Int  Part, 
cm”3  s 

J 

Int  Part  Corr, 
cm”3  s 

K 

AvgCorr, 

am 

L 

CCL/Part, 
!0”12  2 

990412 

Atlas 

84641 

84693 

201.5 

79.0 

2,84 

418 

12.0 

6,14 

6348 

7615 

1,20 

8.06 

84973 

85017 

200.0 

70.0 

2.54 

428 

5,9 

1.74 

2700 

2895 

1.07 

6M 

990415 

Delta 

67017 

67023 

210.6 

65.0 

2.24 

460 

27.0 

1.36 

1684 

2896 

1.72 

4.70 

67936 

67944 

212.1 

191.0 

6.53 

341 

12.0 

3.08 

2558 

4261 

1.67 

7.23 

67947 

67952 

213.5 

188.0 

6.38 

344 

9.0 

1.39 

1243 

1810 

1.46 

7.68 

68135 

68148 

214,0 

203.0 

6.82 

338 

9.0 

2.61 

2348 

3763 

1,60 

6.94 

68145 

68149 

213.8 

202.0 

6.85 

338 

2.3 

0.28 

371 

428 

1.15 

6,54 

68379 

68387 

214.5 

199.0 

6.73 

340 

8,8 

2.09 

2015 

3091 

1.53 

6.76 

68389 

68395 

214.2 

204.0 

6.90 

338 

2.4 

0.32 

458 

505 

1.10 

6.33 

68408 

68411 

213.8 

206,0 

6.99 

336 

1.6 

0.14 

185 

188 

1.02 

7,43 

68523 

68535 

218.2 

225.0 

7.48 

334 

2,5 

0.50 

612 

688 

1.12 

7.27 

68536 

68549 

217.5 

221.0 

7,37 

335 

4.0 

1.04 

1324 

1430 

1.08 

7.27 

68550 

68562 

216.5 

214.0 

7,16 

336 

5.7 

2,02 

2210 

2710 

1,23 

7.45 

68742 

68750 

212.1 

194.0 

6.63 

339 

7.6 

2.06 

1920 

2635 

1.37 

7.82 

68752 

68765 

212,1 

198.0 

6,73 

337 

5,7 

1,27 

1396 

1717 

1.23 

7.40 

68767 

68787 

213.1 

204.0 

6.93 

336 

6.7 

2.00 

2374 

3045 

1.28 

6.57 

68956 

68974 

214.0 

206.5 

7.03 

336 

5.3 

2.02 

2385 

2962 

1.24 

6,82 

69146 

69158 

212.5 

193.5 

7.03 

335 

8.5 

2.08 

2306 

2942 

1.28 

7.07 

69166 

69206 

212.6 

195.5 

6.85 

337 

4.5 

3.57 

4155 

4741 

1.14 

7.53 

69357 

69369 

213.6 

205.0 

6.95 

337 

6,1 

1.45 

1720 

2234 

1.30 

6,49 

69591 

69601 

212.4 

192.0 

6.53 

341 

7.6 

1.90 

1841 

2376 

1.29 

8.00 

69619 

69652 

212.7 

201.0 

6.83 

335 

6.8 

4.28 

5193 

6352 

1.22 

6,74 

69807 

69822 

213.3 

201.2 

6.84 

335 

4,5 

1.35 

1731 

2104 

1.22 

6,42 

69855 

69860 

212.6 

192.5 

6.58 

340 

3,5 

0.53 

602 

683 

1.13 

7 16 

990924 

Athena 

66376 

66384 

208.1 

70.0 

2,44 

445 

8,5 

0.40 

1052 

1610 

1,53 

2.45 

66728 

66737 

207.8 

68.8 

2.40 

447 

4.7 

0.33 

1916 

2480 

1,29 

1.31 

67073 

67080 

207.2 

66.5 

2.33 

450 

3.8 

0.20 

1204 

1627 

1.35 

1.22 

67363 

67382 

207.0 

68.6 

2.40 

445 

5,3 

0.44 

1975 

2629 

1.33 

1,69 

67721 

67746 

206.9 

65.9 

2.30 

451 

4.3 

0.53 

2846 

3759 

1.32 

1.40 

68311 

68324 

201.5 

93.2 

3.35 

397 

2.5 

0.31 

1219 

1394 

1,14 

2.24 
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ppm  at  a  rate  of  ~0.01  STP  L  min*1.  Ambient  air  is 
compressed  to  ~500  hPa  with  a  modified  one-stage  Tef¬ 
lon-valve  diaphragm  pump  (KNF  Neuberger  N-026),  and  a 
portion  of  this  flow  is  passed  through  the  second  cell  (the 
“sample”)  at  a  rate  of  ~0.4  STP  L  min”1  providing  a 
refresh  rate  of  about  Is”1. 

[12]  The  absolute  pressure  in  the  reference  cell  and  the 
pressure  differential  between  the  two  cells  are  controlled 
with  custom-built  electronics  modeled  after  a  commercial 
controller  (Model  1250,  MK3  Instruments,  Andover,  Mass.). 
The  detector  is  packaged  in  a  sealed  vessel,  into  which  the 
flow  through  the  reference  cell  is  drained  such  that  the  entire 
contents  of  the  vessel  are  held  at  the  same  pressure.  Over  the 
typical  duration  of  a  plume  encounter  (tens  of  seconds) 
variations  of  temperature  and  pressure  in  the  box  and  the 
two  detector  cells  are  negligible  relative  to  electronic  noise. 
Consequently,  atmospheric  variability  and  electronic  noise 
dominate  the  variations  in  signal  reported  here. 

[13]  For  simplicity,  and  to  avoid  undesirable  shifts  during 
plume  encounters,  the  temperature  of  the  instrument  was  not 
controlled;  rather,  it  was  allowed  to  reach  a  self-heated 
value  of  about  30°C,  Extensive  calibrations  in  flights  not 
shown  here  demonstrated  that  the  instrument  response 
varies  by  less  than  1%  at  a  constant  pressure,  and  is  well 
described  by  the  ideal  gas  law.  For  the  data  reported  here  the 
precision  of  the  instrument  was  better  than  1  ppm  at  >1  Hz 
response  and  was  improved  with  each  flight  by  attenuating 
various  sources  of  electronic  noise. 

[14]  The  accuracy  of  the  C02  measurements  depends  on 
that  of  the  standards  used  to  calibrate  the  instrument.  For 
these  flights,  with  emphasis  on  fast  measurements  of 
enhancements  in  narrow  plumes,  the  standards  were  only 
characterized  to  a  few  ppm,  such  that  absolute  C02  abun¬ 
dances  are  accurate  only  to  this  level.  Enhancements  of  C02 
observed  during  plume  encounters  are  accurate  to  the  level 
of  uncertainty  of  the  response  (or  “span”)  of  the  detector, 
which  we  estimate  to  be  better  than  ±5%  at  2a.  Mixing  ratios 
of  ambient  C02  varied  between  365  and  370  ppm  during 
these  flights;  In  this  paper  we  are  concerned  only  with  plume 
enhancements  In  C02  above  ambient  values.  The  accuracy 
of  those  enhancements  is  limited  by  the  knowledge  of  C02 
in  the  air  around  the  plume  wake,  which  is  estimated  from 
observations  of  C02  immediately  before  and  after  plume 
ingress  and  egress,  respectively.  Except  where  noted,  the 
overall  2cr  uncertainty  of  the  plume  enhancements  of  C02  is 
the  larger  of  ±0.5  ppm  or  ±5%. 

[is]  Number  and  size  distributions  of  particles  with 
diameters  in  the  range  0.37  to  4.0  pm  were  measured  by 
laser  particle  counters  (LPCs)  (Met  One  Model  237B) 
repackaged  for  the  WB-57F  [Ross  et  al. ,  1999].  This  size 
range  likely  includes  most  of  the  mass  of  alumina  emitted 
by  the  rocket,  although  there  is  no  consensus  on  the  relative 
mass  in  the  submicron  portion  of  the  size  distribution,  as 
noted  earlier.  Soot  particles  produced  by  the  kerosene  fueled 
engines  of  the  ATLAS  IIAS  and  Delta  II  rockets  are  likely 
too  small  (~0, 02— 0,04  pm)  to  be  detected.  The  sizes  of 
particles  measured  in  flight  are  determined  relative  to 
laboratory  calibrations  using  polystyrene  latex  spheres. 
The  residence  time  of  sample  air  between  the  inlet  and 
the  LPC  detector  volume  is  about  1 .5  s. 

[16]  For  the  Atlas  IIAS  and  Delta  II  plume  encounters,  a 
single  LPC  measured  the  total  particle  population.  For  the 


Athena  II  plume,  two  LPCs  were  employed  in  series,  one 
to  measure  total  particles  and  the  other  to  measure  par¬ 
ticles  remaining  after  the  airflow  was  heated  to  I98°C. 
Because  water,  HCI,  and  HNO3  evaporate  at  much  lower 
temperatures,  the  heated  LPC  channel  is  believed  to 
represent  SRM  alumina  and  a  small  amount  of  aluminum 
chlorides,  whereas  the  unheated  LPC  channel  represents 
the  total  aerosol  including  volatile  components  that  survive 
the  transit  between  the  instrument  inlet  and  the  first  LPC. 

[17]  Enhancements  of  large  particles  within  rocket  plumes 
are  typically  several  orders  of  magnitude  larger  than  back¬ 
ground  concentrations,  such  that  the  particle  measurements 
are  highly  precise.  For  these  flights,  number  densities  of  large 
particles  were  obtained  every  0.6  s,  whereas  size  distribu¬ 
tions  were  determined  whenever  counts  from  5000  particles 
were  accumulated.  At  high  particle  number  densities,  particle 
counters  are  subject  to  saturation  and  undercounting  (‘coin¬ 
cidence  losses’).  Detailed  laboratory  tests  with  the  LPCs 
indicate  that  the  correction  is  small  (less  than  10%)  at 
densities  below  150  cm”3,  but  exceed  100%  at  particle 
densities  of  400  cm”3  and  larger.  Although  densities  larger 
than  350  cm”3  were  observed  mainly  during  the  earliest 
plume  encounters,  some  correction  was  necessary  for  all  of 
the  data.  For  the  results  reported  here,  this  correction  was 
determined  from  data  obtained  in  the  plume  of  the  Delta  II 
rocket,  as  described  later.  The  most  significant  uncertainties 
in  the  particle  measurements  are  due  to  sampling  effects, 
primarily  inertial  enhancements  around  the  WB-57F  air¬ 
frame  and  instrument  inlet  and  impaction  losses  in  transfer 
lines. 

[is]  Water  vapor  measurements  were  obtained  only  on 
September  24,  1999,  in  the  Athena  II  plume  wake  with  the 
JPL  Laser  Hygrometer  [May,  1998].  The  instrument  features 
a  near-infrared,  tunable  diode  laser  (TDL)  spectrometer  and 
an  open-path  optical  cell  external  to  the  aircraft,  minimizing 
sampling  issues  such  as  wall  losses  or  evaporation  of  volatile 
particles.  The  JPL  Laser  Hygrometer  has  a  precision  of  0.030 
ppm  for  1.4-s-averaged  data,  and  an  accuracy  of  ±5%  at  19 
km  altitude  [Kiev  et  al. ,  2000].  Pressures  and  temperatures 
were  measured  at  the  nose  of  the  WB-57F  using  a  Weston 
7881  T  Series  vibrating  cylinder  with  an  accuracy  of  better 
than  ±1  hPa,  and  a  Rosemount  thermistor  with  an  accuracy  of 
better  than  ±1  K. 

[t 9]  Because  we  are  interested  in  the  point-by-point  cor¬ 
relation  between  C02  and  the  other  measurements,  the  data 
sets  must  be  appropriately  averaged  and  aligned.  The  particle 
measurements  are  inherently  faster  than  those  of  C02,  and  so 
a  five-point  running  mean  has  been  applied  to  the  raw  5  Hz 
C02  data.  To  account  for  the  longer  residence  time  of  air  In 
the  C02  sample  lines  compared  to  the  particle  and  water 
instruments,  the  C02  data  were  shifted  backward  In  time  by 
a  few  seconds  based  on  the  simultaneous  5-Hz  measure¬ 
ments  of  CIO  (D.  Toohey,  unpublished  data,  1999).  The 
results  are  then  merged  back  into  a  spreadsheet  with  a 
timebase  of  0.6  s  as  determined  by  the  particle  measure¬ 
ments.  Thus,  It  Is  important  to  note  that  for  the  particle-C02 
correlations,  two  consecutive  C02  data  points  are  not 
entirely  independent.  For  the  H20-C02  comparisons,  the 
data  rates  are  comparable  and  the  ordered  pairs  can  be 
considered  as  independent  measurements  (that  is,  only  every 
other  C02  data  point  is  included  in  the  analysis  because  the 
water  measurements  are  reported  every  1.4  s). 
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4.  Results 

4.1.  Correlation  of  C€>2  and  Particles 

[20]  Time  series  of  C02  and  particle  number  densities 
observed  during  the  sixth  pass  through  the  plume  of  the 
Delta  II  rocket  are  shown  in  Figure  1.  As  was  typical  for  all 
the  plume  encounters  for  each  of  the  rockets,  the  two 
quantities  are  highly  correlated*  The  C02  signal  does  not 
respond  as  rapidly  to  small-scale  features  as  does  that  from 
the  particle  instrument  (e.g.,  the  peak  at  68770  s),  as 
expected*  Such  a  tight  correlation  indicates  that  the  distinct 
exhaust  plumes  from  the  Delta  II  rocket  (a  single  kerosene 
fueled  engine  surrounded  by  three  SRMs)  are  well  mixed  by 
turbulent  processes  by  the  time  of  the  WB-57F  encounters* 

[21]  A  scatterplot  of  C02  mass  and  particle  number 
densities  for  the  entire  set  of  low-altitude  Delta  II  plume 
observations  is  shown  in  Figure  2.  The  data  are  not  ordered 
with  respect  to  time,  indicating  that  the  relationship 
between  C02  and  particle  number  is  conserved  as  the 
plume  slowly  expands  and  mixes  with  background  air. 
The  scatter  in  Figure  2  can  be  explained  by  small  differ¬ 
ences  in  the  response  times  of  the  instruments  and  the  lower 
precision  of  the  C02  measurements, 

[22]  The  main  feature  of  Figure  2  is  the  nonlinearity  of 
particle  number  densities  with  increasing  C02.  Laboratory 
tests  of  the  LPC  with  sodium  chloride  particles  found 
coincidence  losses  large  enough  to  explain  this  nonlinearity. 
These  tests,  which  are  referenced  to  an  independent  particle 
counter,  are  somewhat  uncertain  at  the  highest  particle 
densities  observed  in  dense  plume  regions.  Therefore,  based 
on  the  correlation  shown  in  Figure  2,  we  derive  a  third-order 
polynomial  function  that  best  reduces  the  nonlinearity, 

Pc  =  p0  -  L0  x  ID*4  Pj  +  5,0  x  10”6  P*f  (4) 

where  P0  and  Pc  are  the  uncorrected  and  corrected  particle 
number  densities.  A  similar  order  is  sufficient  to  explain  the 
curvature  observed  in  the  laboratory  tests.  Although  not 


Figure  1.  Traces  of  C02  mass  density  enhancements 
(solid  line)  and  uncorrected  large  particle  number  densities 
(dashed  line)  in  the  plume  of  the  Delta  II  rocket  at  ~11  km. 
The  high  degree  of  correlation  on  short  timescales  Is  typical 
of  the  measurements  in  other  passes  and  for  other  rocket 
types. 


Figure  2.  Scatterplot  of  enhancements  of  C02  versus 
enhancements  of  large  particles  for  all  of  the  low-altitude 
(-11-12  km)  passes  through  the  Delta  II  plume  (67936- 
69860  s).  Data  from  outside  the  plumes  have  been  removed 
for  clarity* 


shown  here,  application  of  this  empirical  correction  factor  to 
the  Atlas  HAS  and  Athena  II  plume  data  also  markedly 
improves  the  linearity  of  those  correlations.  On  this  basis, 
we  conclude  that  the  correlation  nonlinearity  is  an 
instrumental  effect,  rather  than  a  property  of  the  plumes, 
and  that  both  the  C02  and  particle  instruments  behaved 
consistently  on  all  three  flights.  The  coincidence  loss 
corrections  derived  using  equation  (4)  increase  the  total 
number  of  particles  observed  in  the  Atlas  HAS,  Delta  II,  and 
Athena  II  plumes  by  16,  29,  and  32%,  respectively.  The 
values  of  the  average  corrections  for  each  plume  encounter 
are  summarized  in  Table  2, 

[23]  Table  2  also  presents  the  quantities  of  C02  mass  and 
particle  number  densities  integrated  across  each  plume 
encounter.  These  integrated  quantities  are  useful  because 
they  reduce  random  errors,  eliminate  scatter  due  to  differ¬ 
ences  in  instrument  time  response  and  data  reporting  times, 
and  reduce  the  effects  of  incomplete  mixing  in  the  plumes. 
There  is  a  high  degree  of  consistency  from  encounter  to 
encounter  for  each  rocket,  especially  for  the  Delta  II, 
Figure  3  show's  the  correlation  between  integrated  mass  of 
C02  and  Integrated  number  of  particles  for  all  the  low- 
altitude  plume  encounters  of  the  Delta  II  rocket.  Although 
covering  a  large  range  in  the  magnitudes  of  these  integrated 
quantities  due  to  shearing  and  mixing  processes,  the  results 
using  corrected  particle  data  fall  on  a  line  with  a  correlation 
coefficient  of  0,99,  whereas  results  based  on  the  uncorrected 
particle  data  are  significantly  curved. 

[24]  C02  mass-to-partiele  number  ratios  integrated  across 
each  discreet  plume  segment  are  plotted  in  Figure  4  against 
the  particle  coincidence  loss  corrections  averaged  across 
each  plume.  There  is  a  clear,  systematic  increase  in  these 
ratios  with  increasing  average  coincidence  loss  correction 
for  the  uncorrected  particle  number  data,  whereas  the  values 
using  corrected  particle  number  data  are  randomly  scattered 
about  an  average  of  —7,2  x  10“ 12  g  C02  particle-1.  Further 
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Integrated  particle  number  density  (cm*3  s) 


Figure  3.  Correlation  plot  of  C02  mass  density  enhance- 
ment  versus  particle  number  density  enhancement,  inte¬ 
grated  across  each  discrete  plume  segment  listed  in  Table  2 
for  the  Delta  II  rocket  Open  circles  are  determined  from  the 
“raw”  particle  number  density  measurements,  whereas  the 
solid  circles  represent  those  particle  data  that  are  corrected 
for  coincidence  losses.  The  outlier  at  integrated  particle 
density  ^3000  em“3  is  the  first  encounter  at  19  km. 


inspection  of  Table  2  reveals  several  trends  in  die  data.  First, 
for  the  measurements  at  lowest  densities  (19  km)  there  is  a 
clear  ordering  of  the  C02/partIcIe  number  ratio  of  about  7.5, 
4.7,  and  L6  (1(T12  g  C02  particle*1)  (4.7  to  2.9  to  1)  for 
the  Atlas  IIAS,  Delta  II,  and  Athena  II,  respectively. 
Assuming  well  mixed  plume  wakes  at  19  km,  the  CO 2/ 
alumina  mass  ratios  for  Atlas  IIAS,  Delta  II,  and  Athena  II 
are  expected  to  be  5.9,  2.7,  and  1.1  (5.4  to  2.5  to  1), 
respectively,  similar  to  the  observed  ordering.  This  suggests 
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Figure  4.  Seatterplot  of  the  ratios  of  C02  mass-to-particle 
number  for  the  low-altitude  encounters  with  the  Delta  II 
plume  as  a  function  of  the  average  coincidence  loss 
correction.  Symbols  are  as  in  Figure  3. 


that  the  alumina  particles  emitted  by  the  different  SRMs 
used  by  the  three  rockets  have  similar  physical  character¬ 
istics  (e.g.,  composition,  diameter,  and  density). 

[25]  Second,  the  data  show  that  C02  mass  relative  to  the 
number  of  particles  is  correlated  with  altitude  (or  ambient 
pressure),  whereas  the  ratio  of  the  mass  of  C02  to  mass  of 
alumina  Is  not  expected  to  vary  with  altitude.  With  the 
exception  of  the  first  encounter  of  the  Athena  II  plume,  the 
data  reveal  a  clear  trend  toward  fewer  particles,  relative  to 
C02,  with  decreasing  altitude.  This  trend  is  most  striking  in 
the  Delta  II  data,  where  the  ratio  of  mass  of  C02  to  mass  of 
alumina  is  expected  to  be  greater  at  19  km  than  at  11-12 
km,  whereas  the  data  present  the  opposite  trend.  Thus, 
assuming  the  particles  have  similar  size  distributions  and 
densities,  there  are  consistently  more  particles  at  high 
altitudes  compared  to  low  altitudes.  Because  we  expect 
C02  to  be  a  conservative  tracer  of  propellant  consumption, 
this  trend  must  be  the  result  of  altitude-related  variations  In 
particle  production  mechanisms  or  in  the  detection  effi¬ 
ciency  of  the  LPCs.  This  conclusion  is  supported  by 
detailed  analysis  of  the  Delta  II  plume  particle  size  distri¬ 
butions  presented  below. 

4.2.  Particle  Size  Distributions 

[26]  The  size  distributions  measured  In  each  of  the 
plumes  are  broadly  consistent.  Figure  5  shows  the  average 
of  the  differential  number  densities,  dN/dD,  versus  particle 
diameter  for  each  of  the  three  rockets  (after  discarding  the 
first  encounter  to  minimize  bias  due  to  the  relatively  large 
coincidence  loss  corrections).  To  highlight  differences  in  the 
shape  of  the  distributions,  we  have  normalized  the  results  to 
the  values  at  0.91  pm.  There  is  a  mode  of  particles  near  1,1 
pm  and,  for  the  total  LPC  channel  data,  an  increase  in 


Figure  5.  Average  large-particle  size  distributions  ob¬ 
served  for  each  of  the  three  rocket  types  (excluding  the  first 
pass  through  each  of  the  three  rockets  where  coincidence 
loss  corrections  were  largest),  normalized  to  the  value  at 
0.91  pm.  The  total  particle  distributions  from  the  three 
rockets  are  virtually  indistinguishable,  whereas  the  non¬ 
volatile  channel  from  the  Athena  II  plume  intercepts  reveal 
a  shift  to  smaller  diameters  and  a  relative  “loss”  of  0,37  pm 
particles. 
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differential  number  density  at  the  smallest  particle  diameters 
detected  by  the  LPC,  suggestive  of  an  additional  mode 
outside  of  the  LPC  window.  There  is  a  hint  of  another,  more 
variable,  mode  near  3  pm  for  all  of  the  rockets,  as  reported 
in  a  previous  study  [Cofer  et  al,  1987].  Otherwise,  the 
distributions  look  remarkably  similar. 

[27]  Focusing  on  the  most  extensive  data  set  with  the 
smallest  potential  sampling  errors,  the  Delta  II  low  altitude 
data,  we  determine  the  volume-averaged  diameter  from 
measured  size  distributions  using  the  equation 


D(cm)  = 


te^N(r)v(r)6(r) 
^EN(r)s(r)  ’ 


(5) 


where  v(r)  is  the  volume  of  a  sphere  and  N(r)  is  the  size 
distribution  measured  in  size  bins  of  width  6(r).  Gravita¬ 
tional  settling  and  coagulation  of  particles  smaller  than  10  pm 
diameter  are  slow  on  the  timescales  of  these  measurements 
such  that  any  changes  In  particle  sizes  must  be  due  to 
condensation,  evaporation,  or  evolution  of  particle  density. 
A  volume-weighted  mean  particle  diameter  of  1.5  pm  is 
obtained  If  this  average  size  distribution  is  extrapolated  out 
to  6  pm  assuming  a  lognormal  distribution.  This  result  is 
similar  to  those  obtained  In  the  plumes  of  large  SRMs  in  the 
midtroposphere  [Cofer  etal,  1987]  and  stratosphere  [Ross  et 
ah ,  1999].  Based  on  the  C02/particle  ratios  in  Table  2  and 
the  known  CCh-to-alumina  mass  ratio  of  1.8  for  the  Delta  II 
emission  at  11-12  km,  the  particles  would  have  a  mean 
density  of  2  g  cm  3. 

[2s]  There  is  considerable  uncertainty  regarding  the  den¬ 
sity  of  alumina  particles  emitted  by  SRMs  in  the  strato¬ 
sphere  [e.g,.  Betting,  1997].  The  value  we  derive  is  at  the 
lower  end  of  the  range  assumed  by  Turco  et  al  [1982]  and 
reported  by  Strand  et  al,  [1981],  but  nearly  identical  to  the 
most  likely  value  deduced  by  Cofer  et  al  [1987].  Because 
inertial  enhancement  effects  wrou!d  shift  the  apparent  size 
distributions  to  larger  sizes,  and  smaller  particles  may  be 
composed  of  volatile  materials  (see  below),  the  density 
determined  from  our  measurements  Is  likely  a  lower  limit 
to  the  true  value.  Even  so,  we  conclude  that  the  particles  are 
composed  primarily  of  alumina  and  that  the  single  1.1  pm 
particle  mode  can  account  for  all  of  the  alumina  emitted  by 
the  SRMs.  Given  the  uncertainty  in  the  effective  alumina 
density,  these  data  can  not  completely  exclude  some  mass 
being  In  modes  smaller  than  0.3  pm.  Because  the  density  of 
alumina  likely  is  not  much  smaller  than  2  g  cm-3,  we 
conclude  that  submicron  particles  probably  do  not  contain 
more  than  a  few  tens-of-percents  of  the  total  mass  of 
alumina  emitted,  eliminating  the  most  extreme  value  that 
has  been  reported  previously  [Dao  et  al,  1997], 

[29]  Assuming  that  the  density  of  the  particles  observed  at 
19  km  Is  also  ~2  g  cm”3,  the  C02/particle  ratios  in  Table  2 
and  the  known  emission  indices  imply  a  volume-weighted 
mean  particle  diameter  of  1. 1  pm  for  all  three  rockets.  This 
value  Is  somewhat  smaller  than  that  of  1 .5  pm  based  on  the 
size  distributions  shown  in  Figure  5.  This  difference  can  be 
explained  In  several  ways.  The  mean  density  of  alumina 
particles  could  be  60%  smaller  at  19  km  than  at  1 1  - 12  km, 
possibly  because  particle  density  varies  inversely  with  size 
[Belting,  1997]  or  is  influenced  by  the  presence  of  low 
density  compounds.  Alternatively,  there  could  be  a  signifi¬ 
cant  inertial  enhancement  of  particles  larger  than  1  pm  that 


Figure  6.  Scatterplot  of  observed  total  water  vapor 
volume  mixing  ratios  (In  parts  per  million)  versus  C02 
volume  mixing  ratio  enhancements  (in  parts  per  million)  for 
the  first  five  passes  through  the  plume  of  the  Athena  II 
rocket  (66376-67746  s).  Data  outside  of  the  plumes  have 
been  removed  for  clarity. 


increases  with  lower  pressure  and  larger  WB-57F  airspeeds. 
Finally,  the  relative  number  of  the  smallest  particles  detect¬ 
able  by  the  LPC  may  be  greater  at  19  km  than  at  lower 
altitudes.  We  discuss  some  important  implications  of  these 
possibilities  in  the  next  section. 


5.  Discussion 

[30]  To  first  order,  the  particles  emitted  by  the  Atlas  HAS, 
Delta  II,  and  Athena  II  rockets  are  characterized  by  a  single 
mode  of  1 .1  pm  with  a  density  of  about  2  g  cm”3,  consistent 
with  previous  observations  of  alumina.  In  addition,  at  19  km 
the  comparative  number  of  particles  relative  to  C02  for  each 
rocket  is  consistent  with  the  trend  In  comparative  emissions 
of  alumina  relative  to  C02  with  rocket  type  (Table  1).  Other 
aspects  of  the  data  require  considerable  elaboration  to  this 
overall  picture,  however.  We  showed  above  that  the  number 
of  particles  relative  to  C02  (i.e.,  the  particle  number 
emission  Index)  tends  to  increase  with  decreasing  ambient 
air  density.  It  is  reasonable  to  expect  this  to  be  due,  at  least 
in  part,  to  nonrepresentative  sampling  of  larger  particles  as  a 
result  of  inertial  enhancement  effects  that  increase  with 
altitude  (e.g.,  lower  ambient  air  densities  and  greater  WB- 
57F  airspeeds).  Detailed  characterization  of  the  airflow 
around  the  WB-57F  wing  pod  (where  the  LPC  was  located) 
and  the  airflow  in  the  LPC  inlet  would  be  required  In  order 
to  evaluate  the  magnitude  of  the  inertial  enhancements, 
which  is  beyond  the  scope  of  this  study.  Thus,  we  can  only 
assume  that  some  part  of  the  data  is  explained  by  such  an 
effect. 

[31]  There  is  a  more  intriguing  explanation  of  the  altitude 
variations  of  the  particle  characteristics.  Volatile  compounds 
condensed  onto  the  alumina  particles  would  modify  their 
size  and  optical  properties  and,  consequently,  their  ability  to 
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be  detected  by  the  LPC.  A  variety  of  volatile  compounds, 
such  as  H20,  HC1,  and  HN03,  are  known  constituents  of 
SRM  plumes.  Based  on  measured  C02  and  the  expected 
emission  ratios,  H20  and  HC1  enhancements  In  the  Athena 
II  plume  would  be  about  5-10  ppm  and  1-2  ppm, 
respectively,  and  the  relative  humidity  in  some  regions  of 
the  exhaust  plume  would  likely  exceed  100%,  HN03 
enhancements  of  20-40  ppb  were  measured  by  Neuman 
ei  ah  [2000]  In  the  Athena  II  plume  wake.  With  alumina 
serving  as  condensation  nuclei,  the  conditions  for  volatile 
condensation  are  expected  in  the  Athena  II  plume.  Direct 
evidence  for  condensed  volatile  compounds  is  provided  by 
measurements  of  water  vapor  mixing  ratios  in  the  Athena  II 
plume,  which  are  shown  in  Figure  6  as  a  function  of  C02 
mixing  ratios.  Although  highly  correlated  with  C02,  water 
vapor  abundances  exhibit  a  shortfall  of  about  0,5  to  1.5  ppm 
over  the  range  of  observed  C02  enhancements  compared  to 
the  expected  H20  to  C02  molar  emission  ratio  of  1.9, 
Condensation  of  this  quantity  of  H20  onto  micron-sized 
alumina  particles  at  typical  plume  number  densities  (~102 
cm  3)  would  increase  the  apparent  particulate  volume  by 
about  one  third.  Alternatively,  the  same  quantity  of  H20 
condensed  onto  0.3  pm  alumina  particles  could  activate 
more  than  500  cm“3  so  that  they  would  become  detectable 
by  the  LPC,  assuming  they  did  not  completely  evaporate 
during  the  ~1.5  s  transit  from  the  inlet  to  the  LPC  detection 
volume. 

[32]  A  comparison  of  the  Athena  II  plume  nonvolatile  and 
total  LPC  size  distributions  in  Figure  5  shows  that  the 
relative  number  of  the  smallest  particles  was  indeed  signifi¬ 
cantly  reduced  after  the  sample  was  heated,  presumably  due 
to  evaporation  of  all  water  condensed  on  these  particles. 
There  was  also  a  shift  in  apparent  size  toward  larger  particles 
for  the  total  LPC  distribution,  which  cannot  be  explained  by 
condensation  of  only  1  ppm  of  water  vapor.  However, 
because  the  optical  scattering  properties  of  ice-coated  alu¬ 
mina  are  significantly  different  from  those  of  bare  alumina, 
the  LPCs  could  be  overestimating  the  diameters  of  ice- 
containing  particles.  Whether  the  result  of  one  or  a  combi¬ 
nation  of  these  effects,  the  water  vapor  shortfall  and  the 
different  size  distributions  reported  by  the  total  and  non¬ 
volatile  LPC  channels  provide  strong  evidence  of  condensed 
water  in  the  Athena  II  plume. 

[33]  One  difficulty  with  the  conclusion  that  the  Athena  II 
particles  contained  significant  quantities  of  volatile  water  is 
that  pure  ice  is  unstable  at  the  temperatures  and  water  vapor 
partial  pressures  in  these  plumes.  Micron-sized  pure  ice 
particles  produced  in  the  plumes  would  have  evaporated  in 
~102  s,  before  the  WB-57F  plume  encounters  [Haynes  ei 
ah,  1992].  However,  H20  equilibrium  vapor  pressures  over 
admixtures  of  ice  and  nitric  acid  are  significantly  reduced 
[Hanson  and  Mauersberger ;  1988]  at  the  nitric  acid  partial 
pressures  characteristic  of  the  Athena  II  plume  [Neuman  et 
ah ,  2000],  It  is  interesting  to  speculate  how  only  0,05  ppm 
of  total  HNO3  could  inhibit  evaporation  of  nearly  1  ppm  of 
condensed  H20.  An  Intriguing  possibility  is  that  the  HN03 
has  combined  with  H20  to  form  a  few  monolayers  (e.g,, 
~10t5  molecules  cm“2)  of  a  liquid  or  solid  (e.g.,  nitric-acid 
trihydrate,  or  HN03*3H20).  Such  a  coating  on  several 
hundred  submicron  particles  would  consume  only  a  few 
ppb  of  gas-phase  HN03,  As  argued  previously  [Deshler  et 
aly  1994,  Peter  ei  ah ,  1994;  Livingston  and  George ,  1999; 


Warshawsky  et  ahy  1999],  such  a  coating  could  significantly 
slow  the  evaporation  of  water  trapped  under  the  coating, 
possibly  even  preventing  sublimation  under  conditions  of 
ice  subsaturation.  Several  modeling  studies  have  found  that 
die  lifetimes  of  aircraft  contrails  may  be  extended  similarly 
by  uptake  of  HN03  into  liquid  aerosols  (which  may  sub¬ 
sequently  freeze)  [Gleiismann  and  Zellner ;  1999]  or  onto 
surfaces  of  ice  particles  to  form  NAT  [Karcher,  1996], 

[34]  Interestingly,  based  on  the  model  of  Carslaw  et  ah 
[1995]  the  evaporation  points  for  solid  HN03*3H20  (NAT) 
calculated  using  maximum  water  vapor  mixing  ratios  deter¬ 
mined  from  the  C02  measurements  and  expected  emission 
ratios  (Table  1)  were  ^210  and  ~213  K  for  the  Delta  II 
plume  at  19  and  11-12  km,  respectively,  ~205  K  for  the 
Athena  II  plume  (after  the  first  pass),  and  ~207  K  for  the 
Atlas  HAS  plume,  assuming  that  maximum  mixing  ratios  of 
HN03  were  about  50  ppb.  Ambient  temperatures  measured 
during  the  Atlas  HAS  plume  encounters,  the  majority  of  the 
Delta  II  low-altitude  plume  encounters,  and  the  last  encoun¬ 
ter  of  the  Athena  II  plume  was  at  or  below  these  threshold 
values,  such  that  NAT  (or  a  liquid  HN03-H20  mixture) 
were  very  likely  present  during  these  encounters.  Even  more 
interesting,  the  ambient  temperatures  measured  during  the 
19  km  passes  through  the  Delta  II  and  Athena  II  plumes 
were  a  few  degrees  above  NAT  threshold  values.  These 
were  the  intercepts  for  which  the  plumes  were  difficult  to 
locate  visually  by  the  WB-57F  aircrew.  We  note  that 
submicron  HN03*H20  condensates  would  scatter  visible 
light  more  efficiently  than  micron-sized  alumina  particles. 
Therefore,  we  propose  that  plume  visibility  may  be  linked 
to  the  presence  of  condensed  mixtures  of  nitric  acid  and 
water,  and  that  the  persistence  of  condensed  water  ice  in  the 
Athena  II  plume  may  be  the  result  of  a  thin  coating  of  solid 
(i.e.,  NAT)  or  liquid  HN03. 

[35]  Some  condensed  HC1  might  be  expected  In  these 
volatile  particles.  The  partitioning  of  chlorine  In  SRM 
plumes  is  time  variable  and  complicated,  as  gas-phase  and 
heterogeneous  reactions  such  as  reactions  (l)-(3)  convert  it 
Into  active  forms  which  react  quickly  with  ambient  ozone. 
Therefore,  without  time-resolved  measurements  of  HC1  it  is 
difficult  to  speculate  on  its  role  in  the  condensed  phase  for 
these  plumes.  However,  the  same  thermodynamic  model 
[Carslaw  et  ah ,  1995]  suggests  that  it  is  likely  that  imme¬ 
diately  after  the  plumes  form  (when  HC1  mixing  ratios  are  a 
few  ppm),  significant  quantities  of  this  species  would  be  In 
the  condensed  phase.  Such  partitioning  would  very  likely 
accelerate  heterogeneous  conversion  of  C10N02  into  HN03 
via  reactions  (1)  and  (2),  thereby  providing  a  rapid  source 
for  the  HN03  that  would  further  stabilize  ice  particles 
formed  in  the  expanding  plume  wake. 

[35]  Additional  evidence  for  the  presence  of  small  volatile 
particles  In  rocket  plumes  Is  provided  by  multiwavelength 
Lidar  observations  of  SRM  plumes  between  20  and  40  km 
[Dentamaro  et  ah ,  1999].  The  dependence  of  the  plume 
backscatter  with  wavelength  was  found  to  evolve  signifi¬ 
cantly  within  the  first  MG'*  s  after  launch.  Whereas  the 
Rayleigh  scattering  dominated  the  Lidar  signal  shortly  after 
launch,  one  hour  later  the  signal  was  dominated  by  scatter¬ 
ing  from  large  particles.  It  was  suggested  that  this  may  be 
due  to  the  presence  of  evaporating  aqueous  HC1  particles, 
but  the  rate  of  the  observed  changes  was  not  consistent  with 
that  explanation  [Dentamaro  ei  ah ,  1999].  Slowly  evaporat- 
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ing  nitric-acid  coated  particles  can  explain  better  the  Lidar 
baekscatter  observations.  It  was  farther  noted  that  similar 
behavior  was  not  observed  at  40  km,  where  nitric  acid 
hydrates  would  not  be  stable.  We  conclude  that  small 
particles  that  contain  condensed  volatile  materials  likely 
play  Important  roles  in  the  plumes  of  most  rockets. 

6.  Summary 

[37]  Data  obtained  in  the  plume  wakes  of  Atlas  HAS, 
Delta  II,  and  Athena  II  rockets  using  the  WB-57F  high- 
altitude  aircraft  shows  that  C02  mass  and  037  to  4.0  pm 
particle  number  densities  are  highly  correlated.  The  corre¬ 
lation  is  used  to  empirically  correct  the  data  for  coincidence 
losses  in  the  scattering  volume  of  the  laser  particle  counter 
and  to  deduce  that  the  C02  and  particle  instruments  per¬ 
formed  consistently  during  the  three  flights.  Constrained  by 
measurements  of  C02  and  the  expected  emission  ratios  for 
the  three  rockets,  particle  size  distributions  demonstrate  that 
the  mode  of  1 ,1  pm  particles  with  a  density  of  2  g  cm-3  can 
account  for  all  of  the  alumina  mass  emitted  by  the  SRMs  on 
these  rockets.  These  data  also  suggest  that  submicron  modes 
likely  do  not  account  for  more  than  a  few  tens  of  percents  of 
the  total  alumina  mass, 

[3s]  Using  data  obtained  at  19  km  under  similar  ambient 
conditions,  we  show  that  the  observed  differences  in  the 
relative  abundances  of  C02  mass  and  particle  number  for 
different  rockets  are  consistent  with  expected  differences  in 
the  rocket  emissions,  assuming  the  particles  are  composed 
only  of  alumina.  Several  difficulties  complicate  the  picture, 
however,  and  are  difficult  to  explain.  First,  for  constant  rocket 
emission,  the  abundance  of  C02  relative  to  the  number  of 
particles  decreases  with  altitude.  Second,  the  particle  size 
distributions  contain  more  volume  than  can  reasonably  be 
accounted  for  by  alumina  alone.  We  propose  that  these  are  the 
result.  In  some  combination,  of  inertial  particle  sampling 
effects  and  the  presence  of  condensed  water  and  nitric  acid  in 
the  plumes. 

[39]  A  correlation  plot  of  water  vapor  and  C02  in  the 
Athena  II  plume  wake  reveals  a  shortfall,  compared  to  the 
expected  value,  of  about  0.5— 1 ,0  ppm  of  water  vapor  under 
conditions  of  subsaturation  with  respect  to  pure  ice,  a 
shortfall  that  persists  for  more  than  20  min.  In  conjunction 
with  previously  reported  measurement  of  significant  quan¬ 
tities  of  HNO3  in  the  Athena  II  plume,  we  surmise  that 
cocondensed  nitric  acid  and  water  would  likely  have  been 
stable  in  the  conditions  of  the  plume,  and  we  note  that  there 
is  more  than  enough  nitric  acid  to  coat  the  surfaces  of  the 
particles  with  many  monolayers.  We  propose  that  such  a 
coating  limits  the  evaporation  of  water  trapped  In  the 
particle  (likely  with  an  alumina  core),  thereby  reducing 
the  rate  at  which  the  particles  evaporate  in  the  atmosphere. 

[40]  If  confirmed  by  future  measurement  and  modeling 
efforts,  the  presence  of  volatile  particles,  admixtures  of 
nitric  acid  and  water  in  particular.  In  rocket  plume  wakes, 
would  carry  several  significant  implications.  First,  such 
particles  are  known  to  rapidly  convert  reservoir  forms  of 
chlorine  into  ozone-destroying  forms,  such  as  occurs  in  the 
polar  lower  stratosphere.  Existing  models  of  regional  plume 
wake  chemistry  have  not  considered  any  of  the  variety  of 
heterogeneous  reactions  that  would  proceed  on  volatile 
surfaces,  and  such  reactions  may  contribute  to  the  observed 


fast  removal  of  ozone  in  rocket  plume  wakes.  Even  more 
troublesome,  the  lifetimes  of  such  particles,  which  have  not 
been  considered  in  model  studies  of  rocket  impacts,  may  be 
quite  long  in  the  winter  polar  stratosphere  where  conditions 
remain  close  to  NAT  saturation  for  several  months.  Con¬ 
sequently,  high-latitude  rocket  launches  Into  the  polar 
vortex  could  have  Important  implications  for  polar  ozone 
loss. 

[41]  Our  results  also  suggest  that  some  combinations  of 
rocket  engines  may  be  more  harmful  to  stratospheric  ozone 
than  others.  Some  launch  vehicles  combine  SRMs  (emitting 
H20  and  alumina  particles)  and  an  engine  using  N204  and 
hydrazine  compounds.  The  NOx  emission  index  for  hydra¬ 
zine  fueled  rocket  engines  is  predicted  to  be  at  least  twenty 
times  that  for  SRMs  in  the  stratosphere  (unpublished 
combustion  models).  In  the  plume  wakes  of  such  vehicles, 
the  additional  NOx  would  rapidly  form  additional  HNO3  so 
that  this  combination  of  rocket  engines  could  be  a  much 
more  prodigious  source  of  NAT  particles  than  either  engine 
alone.  Finally,  our  results  emphasize  the  significant  gaps  In 
knowledge  regarding  rocket  emissions  and  their  Impact  on 
ozone,  particularly  regarding  particles,  liquid  rocket  emis¬ 
sions,  and  the  unique  interactions  of  the  different  emissions 
used  by  the  various  launch  vehicles  around  the  world. 
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Quantifying  uptake  of  HN03  and  H20  by  alumina  particles 
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[i]  The  goal  of  this  study  is  to  quantify  uptake  of  H20  and  HN03  by  and  estimate  their 
residence  time  on  alumina  particles  in  Athena-2  rocket  plumes.  This  study  uses  in  situ 
measurements  made  in  the  lower  stratosphere  with  the  NASA  WB-57F  high-altitude 
aircraft  on  24  September  1999.  Constraining  the  Atmospheric  and  Environmental 
Research,  Inc.  (AER),  plume  model  with  available  measurements,  we  found  that  (1)  H20 
uptake  coefficient  for  alumina  particles  is  larger  than  3  x  10~4,  (2)  HN03  is  produced  via 
C10N02  +  HC1  — +  Cl2  +  HN03  on  alumina  particles  and  resides  on  their  surfaces  for  5— 

52  min,  and  (3)  alumina  particles  in  the  plume  are  covered  by  100-200  monolayers  of 
adsorbed  H20  and  0.1  —  10  monolayers  of  HN03  under  lower  stratospheric  conditions. 

These  values  are  uncertain  by  at  least  a  factor  of  2.  We  speculate  that  the  H20  coverage 
remaining  on  alumina  particles  accelerates  the  C10N02  +  HC1  -►  Cl2  +  HN03  reaction, 
thus  leading  to  a  larger  than  previously  thought  global  ozone  loss  to  solid-fueled 
rocket  emissions,  especially  if  at  least  several  percent  of  emitted  alumina  mass  are  in 
Submicron  particles.  INDEX  TERMS:  0305  Atmospheric  Composition  and  Structure:  Aerosols  and 
particles  (0345,  4801);  0317  Atmospheric  Composition  and  Structure:  Chemical  kinetic  and  photochemical 
properties;  0340  Atmospheric  Composition  and  Structure:  Middle  atmosphere — composition  and  chemistry; 

KEYWORDS:  alumina  particles,  uptake  of  HN03  and  H20,  rocket  emissions 
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1-  Introduction 

[2]  Launches  of  solid-fueled  rockets  cause  a  particular 
concern  primarily  because  of  emissions  of  chlorine  and 
alumina  particles  directly  into  the  stratosphere,  thus  deplet¬ 
ing  the  ozone  layer  [Prather  et  ah ,  1990;  Ko  et  aL,  1994; 
Jackman  et  aL  1998 ;  World  Meteorological  Organization 
(WMO),  1992].  Recent  analysis  of  solid-fueled  rocket 
plumes  shows  almost  complete  depletion  of  ozone  locally 
[e.g.,  Ross  et  aL ,  2000],  however,  global  implications  of 
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such  sharp  local  ozone  reductions  are  predicted  to  be  small 
[Danilin  et  aL ,  2001a],  On  the  other  hand,  the  heteroge¬ 
neous  reaction  (Rl)  on  alumina  particles: 

(RI )  CIONO2  +  HC1  ^  C12  +  HNOj 

with  the  reaction  probability  of  7  =  0.02  [Molina  et  aL , 
1997]  does  not  play  an  important  role  in  local  ozone 
depletion  [Danilin  et  aL ,  2001a].  However,  it  becomes 
important  on  the  global  scale  by  converting  emitted  and, 
more  importantly  background  HCI  into  short-lived  Cl2 
resulting  in  ozone  depletion  that  would  correspond  to  an 
ozone  depletion  potential  (ODP)  of  alumina  particles  larger 
than  0.2  depending  on  the  location  of  emissions,  size 
distribution  of  alumina  particles,  and  the  value  of  (Rl) 
reaction  probability  [Danilin  et  aL ,  2001b].  Since  the  (Rl) 
rate  is  quite  different  for  pure  [Molina  et  aL ,  1997]  and 
covered  by  H20  [Shi  et  aL,  2001]  alumina,  we  try  to 
quantify  the  uptake  of  H20  and  HN03  by  rocket  plume 
alumina  particles  using  relevant  measurements  obtained 
during  the  Atmospheric  Chemistry  of  Combustion  Emis¬ 
sions  Near  the  Tropopause  (ACCENT)  campaign.  This 
mission  sponsored  by  NASA  and  the  US  Air  Force  took 
place  in  1999-2000  and  was  designed  to  investigate  the 
impact  of  rocket  and  aircraft  emissions  near  the  tropopause. 
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2.  Results  and  Discussion 

2.1,  Measurements 

{3]  During  the  ACCENT  campaign,  emissions  of  four 
different  rockets  were  sampled  by  the  NASA  WB-57F 
aircraft.  We  focused  our  study  on  the  24  September  1999 
flight  because  the  NOAA  chemical  ionization  mass  spec¬ 
trometer  (CIMS)  and  JPL  laser  hygrometer  made  simulta¬ 
neous  measurements  of  HN03  and  H20,  respectively,  in  a 
rocket  plume  during  this  flight  only.  An  Athena-2  was 
launched  at  11,21  am  PST  on  24  September  1999  from 
Vandenberg  Air  Force  Base  In  California.  Figure  1  shows 
the  location  of  the  Athena-2  plume  measurements.  The 
WB-57F  aircraft  made  six  intercepts  of  the  rocket  plume 
with  the  first  five  samplings  occurring  near  70  hPa  (~18.6 
km)  and  the  last  intercept  at  97  hPa  (~16  km).  Plume  ages 
ranged  from  3.7  to  36.2  min  (see  Table  1).  The  burning  of 
1  kg  of  Athena-2  solid  rocket  fuel  resulted  in  1 .27  kg  of 
exhaust  with  emission  Indices  (El)  for  the  main  products 
of  382,  300,  362,  217,  and  2  g/kg(fuel)  for  C02,  H20, 
alumina,  chlorine  (as  HCI),  and  NO,  respectively.  In  this 
study  the  following  measurements  are  used:  C02  [Gates  et 
al >  2002],  H20  [May,  1998],  NOv,  HNO3,  [Popp  et  al , 
2002],  and  particles  [Schmid  et  al,  2002;  Baumgardner  et 
al ,  1996].  Perturbations  of  relevant  species  and  A1203 
particles  in  these  plumes  are  given  by  Popp  et  al  [2002] 
and  Schmid  et  al  [2002],  Fahey  et  al  [1995]  showed  that  it 
is  better  to  use  plume-averaged  values  for  a  quantitative 
analysis,  since  different  instruments  have  different  sam¬ 
pling  rate  and  plume  composition  is  Inhomogeneous. 
Below,  we  follow  this  methodology  for  comparison  of 
plume-averaged  values  for  each  interception  against  our 
model  calculations.  Previous  studies  [Popp  et  al,  2002; 
Herman  et  al,  2001;  Gates  et  al,  2002]  revealed  that  nitric 
acid  and  water  vapor  must  have  condensed  on  particles  in 
the  Athena-2  plumes  in  order  to  explain  their  measure¬ 
ments.  This  finding  is  surprising,  since  for  the  reasonable 
plume  dilution  scenarios  and  ambient  temperature  of  208, 1 
and  200.2  K  at  70  and  97  hPa,  respectively,  H20  and 
HNO3  were  not  saturated  with  respect  to  Ice  or  nitric  acid 
trihydrate  (NAT).  Our  paper  extends  their  analysis  by 
proposing  mechanisms  of  the  HN03  and  H20  condensation 
onto  alumina  particles  under  those  conditions  and  discusses 
implications  of  the  results. 

2.2,  Model 

[4]  In  order  to  address  this  issue,  we  apply  the  Atmos¬ 
pheric  and  Environmental  Research,  Inc.  (AER),  plume 
model  [Danilin  et  al,  2001a],  which  was  initialized  by 
the  Aerospace  model  output  (P.  F.  Zittel,  personal  commu¬ 
nication,  2001)  for  the  Athena-2  emissions  along  its  24 
September  1999  trajectory  with  the  fuel  consumption  rate  of 
14  and  15  g  per  vertical  cm  at  18,6  and  16  km,  respectively. 
The  plume  dilution  rate  was  constrained  by  the  C02 
measurements,  since  carbon  dioxide  Is  chemically  inert 
and  Its  emission  index  is  well-known  (EI(€02)  =  382  g/ 
kg(fuel)),  Since  plume  composition  is  spatially  inhomoge¬ 
neous  and  the  distance  from  the  plume  centerline  for  each 
interception  is  unknown,  we  treat  each  plume  event  indi, 
vidually  in  our  model.  The  model  assumes  a  homogeneous 
distribution  of  species  in  the  plume.  Since  we  study  linear 
processes  here  (such  as  uptake  of  HN03  and  H20  which  Is 


Athena -II  plume  samplings  by  WB-57F  on  990924 


Figure  1.  Location  of  the  Athena-2  plume  measurements 
(red  dot)  and  the  WB-57F  aircraft  flight  track  (green  line) 
near  California  coast  on  24  September  1999. 

linearly  proportional  to  A1203  SAD),  it  is  reasonable  to 
compare  model  results  against  plume-averaged  measure¬ 
ments.  Varying  the  initial  plume  cross  section  and  its 
dilution  rate,  we  match  the  plume  width  and  C02  concen¬ 
tration  at  the  time  of  measurements  (see  Figure  2a),  Since 
the  shape  of  the  plume  sampled  by  the  aircraft  is  unknown, 
for  simplicity  we  assume  that  the  plume  is  circular  with 
diameters  known  from  the  measurements.  This  assumption 
may  introduce  additional  uncertainty  In  our  analysis.  Mix¬ 
ing  with  ambient  air  was  parameterized  as  given  by  Danilin 
et  al  [2001a]  and  the  values  of  ambient  species  used  are  the 
WB-57F  measurements  just  outside  the  plumes. 

[5]  Our  model  starts  calculations  5  s  after  exhaust  when 
plume  temperature  approaches  Its  ambient  values  and  stand¬ 
ard  photochemical  kinetics  [Sander  et  al,  2000]  Is  applica¬ 
ble.  In  this  approach,  we  overlook  important  microphysical 
processes  (like  possible  rapid  H20  condensation  onto  and  its 
evaporation  from  alumina  particles  [Gates  et  al,  2002])  in 
the  near-field.  To  address  this  Issue,  model  sensitivity  studies 
were  made  (see  section  2.5). 

[e]  The  alumina  particles  in  the  plumes  are  treated  as 
non-interactive  spheres  in  40  bins  with  diameters  ranging 
from  0.78  nm  to  6,4  pm.  The  shape  of  the  A12Q3  size 
distribution  is  obtained  by  combining  aerosol  instrument 
measurements  by  the  nucleation-mode  aerosol  mass  spec¬ 
trometer  (NMASS)  and  focused  cavity  aerosol  spectrometer 
(FCAS)  instruments  [Schmid  et  al,  2002]  for  d  <  1.1  pm 
and  from  Multiple-Angle  Aerosol  Spectrometer  Probe 
(MASP)  data  [Baumgardner  et  al,  1996]  for  particles  with 
d  >  1.1  pm,  thus  covering  the  range  from  4  nm  to  20 
micron.  Their  initial  distribution  in  our  model  is  normalized 
in  order  to  be  consistent  with  the  measurements  at  the  time 
of  plume  intercepts  (see  Figure  2c).  Omission  of  possible 
coagulation  of  alumina  particles  is  not  important,  since  it 
does  not  affect  micron-size  particles  (which  determine  the 
alumina  surface  area  density  (SAD)  and  thus  are  the  most 
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Table  1.  Age  and  Diameters  of  Plume  Sampled  and  Plume-Averaged  Measured  and  Calculated  Values  of  Water  and  Nitric  Acid* 


Plume 

Age,  min 

D,  km 

H.O, 

ppmv 

HNO-?,  pphv 

El,  g/kg 

SAD,  [inr/cm3 

Measb 

Model 

Con 

a2 

Measc 

Gasd 

Cond 

uj=0 

t7,=1 

1 

3.7 

1.21 

10.26 

11.62 

1.36 

12.5 

25.6 

3.0 

3.5 

103 

4,0 

2 

9.7 

1.95 

7.68 

8.08 

0.40 

6.2 

19.4 

3.0 

2.5 

52 

15 

3.2 

892 

3 

15.4 

1.15 

7.23 

7.97 

0.74 

5.7 

19,7 

4.9 

2.5 

17 

6 

2.7 

831 

4 

20.3 

2.71 

6.90 

7.78 

0.88 

5.3 

12,8 

2.6 

1,5 

22 

5 

2.0 

461 

5 

26.3 

2.59 

7.51 

8.45 

0.94 

3.4 

19.8 

4.3 

2.6 

31 

7 

2.2 

935 

6 

36.2 

1.91 

8.07 

8.76 

0.69 

2.7 

15.5 

3.0 

2.2 

30 

9 

2.0 

526 

“D,  diameter  Meas,  measured  values;  and  Model,  calculated  values.  The  columns  (con)  show  the  amount  of  condensed  H20  and  HN03  together  with  the 
required  uptake  coefficient  tx  (x  10  )  of  H2O  and  residence  time  (Tj,  in  min)  of  HNO3  on  alumina  Sites,  Derived  EI(NO->)  and  measured  AKO3  SAD  are 
also  listed, 
hGa$~phase  only. 
cGas  and  condensed, 

dotj  =  0. 


Important  in  the  context  of  this  study)  during  the  first  hour 
after  exhaust. 

[7]  Since  the  rocket  near-field  model  calculations  of  NO 
am  quite  uncertain  (P.  F.  Zittel,  personal  communication, 
2001),  we  chose  initial  concentration  of  NO  in  a  such  way 
that  model  NO*  is  consistent  with  the  measured  NOv  at  the 
time  of  plume  intercepts  (see  Figure  2b).  Following  this 
approach,  we  obtained  EI(NO)  shown  in  Table  1  and 
expressed  in  g(N02)/kg(fuel).  Our  EI(NO)  values  are  in 
good  agreement  with  those  reported  by  Popp  et  al  [2002], 
who  used  a  different  approach  analyzing  the  plume  aver¬ 
aged  values  of  C02  and  NOv  in  order  to  derive  EI(NO).  The 
rocket  combustion  models  predict  that  the  initial  partition¬ 
ing  of  chlorine  between  CI2  and  HC1  should  be  35%  and 
65%  by  mass,  respectively  [Zittel,  1994].  In  our  calcula¬ 
tions  we  used  ozone  measurements  [Richard  et  al ,  2001]  in 
order  to  determine  this  partitioning  for  each  plume  inter¬ 
cept.  The  derived  Cl2  fraction  varied  between  5  and  15%. 
Below  we  discuss  uncertainties  related  to  the  initial  chlorine 
partitioning. 

2.3.  HN03 

[s]  Nitric  acid  may  appear  on  alumina  particles  via  both 
heterogeneous  reactions  and  kinetic  uptake  after  collisions 
and  leave  via  desorption.  Thus,  the  temporal  evolution  of 
HN03  condensed  onto  particles  ignoring  plume  dilution  is 
described  by 

dHNO\  SAD  (yvCIONOj  +  (Xj  vt  HNOf  )  HNO \ 

dt  4  0) 

Here  SAD  is  the  surface  area  density  of  A1203  particles,  7  = 
0.02  [Molina  et  al ,  1997],  aj  is  the  uptake  coefficient,  v 
and  V]  are  the  mean  thermal  speeds  of  C!ON02  and  HN03 
molecules,  respectively,  t,  is  the  residence  time  of  HN03 
on  alumina  surface,  superscripts  g  and  c  mean  gas  and 
condensed  phases,  respectively.  Since  the  OH  mixing  ratio 
is  about  0.01-0.1  pptv  according  to  our  model  calculations 
minutes  after  exhaust,  the  gas  phase  production  of  HNO3 
via  N02  +  OH  +  M  — ^  HNO3  +  M  is  several  orders  of 
magnitude  smaller  than  that  of  (Rl)  and  is  omitted  In 
equation  (1).  Contributions  of  N2Os  and  BrON02  hydro¬ 
lysis  to  HN03c  are  much  smaller  than  those  from  (Rl)  in  the 
plume  and  are  neglected.  These  facts  show  that  the  (Rl) 
reaction  is  a  principal  source  of  HN03  in  the  Athena-2 
plume. 


[9]  We  assume  that  the  amount  of  HN03  condensed  per 
unit  of  area  should  be  the  same  for  any  particle  size.  This 
assumption  determines  partitioning  of  HNO/  among  alu¬ 
mina  particles.  Equation  (1)  contains  two  unknown  param¬ 
eters,  namely,  ti  and  oq.  While  we  can  not  determine  them 
simultaneously,  we  can  bracket  the  values  of  Tj  for  the 
smallest  (=0)  and  largest  (=1)  values  of  aj.  Figures  2d-2e 
show  the  evolution  of  HN03  in  gas  and  condensed  phases 
in  the  plume  for  a,  =  1. 

[to]  Comparing  model  calculations  and  measurements, 
one  should  keep  In  mind  that:  (1)  the  JPL  laser  hygrometer 
measures  only  gas-phase  H20  and  (2)  the  NOAA  CIMS 
measures  both  the  gas-phase  and  condensed  HN03.  The 
HN03  instrument  samples  through  a  forward-facing  sub- 
isokinetic  inlet,  which  causes  particles  with  diameter  larger 
than  0. 1  pm  to  be  significantly  enhanced  in  the  sample  flow. 
Thus,  any  condensed  HN03  in  the  plume  is  enhanced  by  the 
factor  EF(r)  for  particles  with  radius  r.  Thus,  knowing  the 
amount  of  HNOf  and  HNOf  in  our  model,  the  apparent 
amount  of  HN03w  as  seen  by  the  CIMS  is  determined  by 

HNCff  =  HNOf  +  HNO?  i .  (2) 

E/=i  m 

Here  EF*  and  n,  are  the  CIMS  enhancement  factor  and  the 
number  of  particles  in  the  /th  bin  (shown  in  Figure  2c), 
respectively.  Figure  2f  compares  apparent  values  of  nitric 
acid  according  to  the  measurements  (black  symbols)  and 
model  calculations  (color  lines),  A  good  agreement  between 
modeled  and  measured  HNO/77  values  is  not  surprising, 
since  we  adjusted  t*  for  obtaining  agreement  The  derived 
values  of  Tj  are  shown  In  Table  1,  ranging  from  17  to  52 
min  and  from  5  to  12  min  for  the  eases  without  and  with 
direct  uptake  of  HN03  vapor  by  alumina  particles,  respec¬ 
tively.  The  derived  values  of  tj  for  plume  1  are  an  outlier 
and  discarded  in  our  conclusions. 

2.4.  H20 

[11]  Since  there  is  no  photochemical  production  of  H2Oc 
in  the  plume,  the  evolution  of  H2Oc  is  described  by 

dHjOf  _  SADol2v2H20^  H2(? 
dt  4 

Again,  we  have  two  unknowns  (a2  and  t2)  and  one 
equation  (3),  We  performed  several  sensitivity  studies  by 
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Figure  3.  The  <*2-1/t  relations  forH20  in  Athena-2  plumes. 


deriving  ol2  for  the  r2  values  of  100,  200,  1000  s,  and  oo. 
The  results  are  presented  in  Figure  3.  The  minimum  values 
of  a2  for  infinite  r2  (i.e.  without  H20  desorption  from 
alumina  particles)  are  shown  In  Table  1  and  are  in  the  range 
of  (2.7- 12.5)  x  10“4.  Extrapolating  our  results  up  to  a2  -  1 
using  the  derived  slope  of  da2/d(l/T)  =  2  s,  one  gets  r2  of 
0,5  s.  These  r2  values  are  of  relevance  for  H20  uptake  by 
liquid  water  in  the  stratosphere  [Li  et  al ,  2001],  Table  1 
shows  that  our  calculations  predict  about  0,5  to  1 .0  ppmv  of 
condensed  H20  in  the  plume.  This  is  in  close  agreement 
with  the  observations  by  Gates  et  al  [2002]  and  Herman  et 
al  [2001], 

2.5.  Sensitivity  Studies 

[12]  We  investigate  the  sensitivity  of  the  results  given  in 
Table  1  to  the  uncertainties  of  (1)  particle  measurements,  (2) 
initial  partitioning  of  chlorine  species  between  Cl2  and  HC1, 
and  (3)  initial  condensed  H2Q  on  alumina  particles.  Table  2 
summarizes  our  findings  assuming  oq  =  1  and  r2  =  1000  s 
for  all  calculations  here. 

[b]  Alumina  particle  measurements  in  the  plume  by 
NMAS3,  FCASS,  and  MASP  may  have  a  total  uncertainty 
of  100%,  which  affects  the  Al2€)3  SAD  and  consequently 
the  derived  values  of  a  and  r  for  H20  and  HNO3.  We  made 
sensitivity  calculations  with  increased  and  decreased  alu¬ 
mina  SAD  by  a  factor  of  2  and  show  our  results  for  doubled 
alumina  SAD  in  Table  2  (4th  and  5th  columns).  The  change 
in  alumina  SAD  leads  to  a  change  of  the  (Rl)  reaction  rate 
(which  is  a  main  source  of  HNOf  in  the  plume)  and  affects 
the  partitioning  of  H20  and  HNO3  between  condensed  and 
gas  phases.  For  larger  A1203  SAD,  larger  values  of  H2Oc 
and  HNOf  and  smaller  values  of  rx  and  r2  are  required  in 
order  to  explain  the  measurements.  In  the  case  of  H20,  an 
increase  (decrease)  of  SAD  by  a  factor  of  two  decreases 
(increases)  the  values  of  t2  by  the  same  factor  (compare  3rd 
and  5th  columns  of  Table  2).  In  the  case  of  HNO3,  its 


Table  2.  Sensitivity  of  T!  and  a2  to  Doubling  A1203  SAD, 
Doubling  Initial  Cl2,  and  Initial  0.5  ppmv  of  Condensed  H20, 
Assuming  =  1  and  t2  =  1000  sa 


Plume 

Baseline 

ri  ot2 

2SAD 

Tt  ot2 

2CLb  r, 

0.5  ppmvc  a? 

1 

103 

14.7 

11.3 

7.3 

13.7 

1.3 

2 

15.0 

9.4 

8.5 

4.7 

4.5 

2.9 

3 

6.3 

11.2 

3.8 

5.6 

3.4 

8,2 

4 

5.4 

12.6 

3.0 

6.3 

2.0 

10.5 

5 

7.2 

10.4 

4.1 

5.2 

2.3 

9.3 

6 

9.1 

10.6 

4.5 

5.3 

6.8 

9.7 

^Sensitivity  of  Tj  is  given  in  minutes,  and  sensitivity  of  values  are 
x  1 0-4. 

bot2  is  the  same  as  in  the  baseline  case. 
ct{  is  the  same  as  in  the  baseline  case. 


photochemical  production  is  involved  which  complicates 
the  dependence  of  Tj  on  SAD.  However,  Tt  is  also  reduced 
by  almost  a  factor  of  2  for  doubled  alumina  SAD  (compare 
2nd  and  4th  columns  of  Table  2).  Similar  conclusions  are 
obtained  from  the  model  runs  with  alumina  SAD  reduced 
by  a  factor  of  2  (not  shown  in  Table  2). 

[14]  Near-field  rocket  combustion  models  [ Zittel ,  1994] 
predict  that  35%  by  mass  of  emitted  chlorine  should  be 
present  as  Cl2  at  18  km  in  the  Athena-2  plume.  Since  our 
plume  calculations,  constrained  by  the  WB-57F  measure¬ 
ments,  show  that  only  5-15%  of  chlorine  is  converted  to 
Cl2  In  order  to  be  consistent  with  the  ozone  measurements, 
we  investigate  how  doubling  of  our  initial  CI2  affects  the 
values  of  Tj  (the  values  of  a2  remain  unchanged).  Initial  Cl2 
and  t j  are  linked,  since  higher  values  of  Cl2  lead  to  larger 
values  of  CIO  and  C10N02,  thus  accelerating  the  formation 
of  HNO3  and  shortening  T|.  The  6th  column  of  Table  2 
shows  that  the  link  between  initial  partitioning  of  the 
emitted  chlorine  species  and  required  residence  time  of 
HNO3  on  alumina  particles  is  quite  pronounced.  However, 
one  should  keep  in  mind  that  this  sensitivity  test  was  a  bit 
hypothetical,  since  the  ozone  values  obtained  in  these 
calculations  were  considerably  smaller  than  those  measured 
by  the  WB-57F.  Nevertheless,  it  is  of  interest  to  perform  this 
sensitivity  test  since  the  discrepancy  between  our  derived 
Cl2  fraction  and  the  fraction  predicted  by  the  rocket  com¬ 
bustion  models  is  larger  than  the  expected  uncertainties  in 
the  rocket  combustion  models.  The  source  of  this  discrep¬ 
ancy  is  unclear  but  our  results  suggest  that  the  combustion 
models  may  significantly  overpredict  the  production  of  Cl2. 
In  order  to  validate  near-field  and  our  far-field  calculations 
of  chlorine  partitioning  in  the  SRM  exhaust,  one  should 
have  simultaneous  measurements  of  CIO  and  HC1.  In  situ 
HC1  measurements  may  be  also  valuable  for  the  study  of 
HC1  interactions  with  alumina  particles  in  rocket  plumes 
and  with  other  aerosol  in  the  background  atmosphere, 

[15]  Near-field  calculations  also  show  a  possibility  for 
H20  condensation  onto  alumina  particles  during  the  first 
few  seconds  after  exhaust.  In  order  to  address  this  possi- 


Flgure  2.  (opposite)  Temporal  evolution  of  plume  diameter  (Figure  2a),  NQ>?  (Figure  2b),  HN03  (Figures  2d-2f,  ai  =  1), 
and  H2Q  (Figure  2g,  t2  =  1000  s)  in  the  Athena-2  plumes.  Measured  size  distribution  of  alumina  particles  and  CIMS 
enhancement  factor  are  shown  in  Figure  2c.  Plume  events  are  color-coded  as  shown  and  measurements  are  shown  by  the 
black  symbols. 
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Figure  4,  Number  of  H20  and  HN03  monolayers  on 
Athena-2  alumina  particles  during  first  two  days  after 
exhaust  Plume  events  are  color  coded  as  shown. 


bility  in  our  model  calculations,  which  overlook  plume 
processes  during  the  first  5  s  downstream,  we  performed 
additional  model  runs  with  the  initial  condensed  phase  H20 
equals  to  0,5  ppmv.  Comparison  of  the  3rd  and  7th  columns 
of  Table  2  shows  the  effect  of  initial  H2Oc  on  a2,  which  is 
larger  for  younger  plumes  and  smaller  for  older  plumes, 
hence  showing  that  the  initial  conditions  become  less 
important  with  time.  For  example,  for  plume  6  the  effect 
of  initial  condensed  H20  leads  to  reduction  of  a2  by  only 
10%.  Since  condensed  water  does  not  affect  the  reaction 
probability  of  (Rl)  in  our  calculations,  the  values  of  Ti 
remain  unchanged  in  this  sensitivity  test 

2.6,  Implications 

[i6]  What  will  happen  with  coverage  of  alumina  particles 
under  background  conditions?  To  answer  this  question,  we 
present  2-day-long  model  runs.  Assuming  a2  =  10” 15  and 
0|  =  2  x  10” 15  cm2  for  the  areas  of  one  molecule  of  H20 
and  HNO3  [Grassian,  2002],  respectively,  and  knowing  the 
evolution  of  A1202  SAD  with  time  from  our  model  calcu¬ 
lations,  we  obtain  the  number  of  H20  (n2  )  and  HN03  (nf  ) 
monolayers  on  alumina  particles  shown  in  Figure  4,  Our 
results  indicate  that  hundreds  of  H20  and  up  to  several 
HN03  monolayers  may  cover  alumina  particles  in  the 
Athena-2  plumes  days  after  emission.  These  asymptotic 
values  show  the  coverage  of  alumina  particles  under  back¬ 
ground  conditions  and  may  be  confirmed  by  the  following 
calculations.  Assuming  that  dH2Cffdi  =  0  under  background 
conditions  and  using  equation  (3),  one  gets 

m  OlHiCF  a2r2Ct2V2H2Og 

',j=^d~  = - 4 - ’  W 

which  gives  n J  =  148  monolayers  for  t2~  1 03  s,  ot2  =  10”3, 
v2  =  49400  cm/s,  and  H20  =  5  ppmv  =  1,2  x  1013  #/ cm3  at 
70  hPa.  Using  the  same  approach  for  HNO3,  one  gets 

w  _  Oi  T|  (yvClONOj  +  €q  VyHNC/j  ) 


The  values  of  nf  range  from  0.3  (o q  =  0)  to  13  (at  =  1) 
monolayers  of  HN03  for  typical  values  of  C10N02  (=0,5 
ppbv)  and  HN03  (=1,4  ppbv)  and  assuming  tj  of  1000  s 
and  300  s  for  at  =  0  and  a2  =  1,  respectively.  These  values  of 
nf  are  consistent  with  those  shown  in  Figure  4.  Assuming 
that  surface  properties  of  alumina  particles  are  the  same  for 
all  solid-fueled  rockets,  our  results  for  nf  and  nf  are  valid 
for  the  Space  Shuttle  too,  which  is  the  main  source  of 
alumina  in  the  global  stratosphere. 

[17]  One  also  can  obtain  from  equations  (1)  and  (3)  the 
amount  of  condensed  H2Oc  and  HNOf,  which  is  of  the  order 
of  105-106  #/cm3  and  102- 103  #/cm3,  respectively,  assum¬ 
ing  the  A1203  SAD  of  4  x  I0~4  pm2/cm3  [Danilin  et  ah , 
2001b,  Figure  5]  at  18  km  at40°N, 

2.7.  Discussion 

[is]  The  hydrophilic  properties  of  the  Athena-2  exhaust 
show  that  there  are  three  stages  In  the  evolution  of  alumina 
particles.  First,  the  fresh  alumina  particles  become  coveted 
very  quickly  primarily  by  H20  during  the  first  minute  after 
emission.  Initial  uptake  coefficients  of  H20  and  HN03  are  of 
the  order  of  10”4-10-3  [Grassian,  2002]  and  the  reaction 
probability  for  (Rl)  is  equal  to  0.02  [Molina  et  aL,  1997]. 
Second,  alumina  particles  are  covered  by  several  hundred 
monolayers  of  H20  containing  small  amounts  of  HNO3  and, 
likely,  other  soluble  acids  (H2S04,  HC1)  [Van  Doren  et  aL , 
1990;  Robinson  et  aL,  1998].  Within  minutes  after  exhaust 
the  alumina  particles  should  behave  like  mildly  acidic  wTater 
droplets  with  a2  and  y  of  an  order  of  0. 1  - 1  [Li  et  aL,  2001 ; 
Shi  et  aL,  2001].  Third,  hours  after  emission,  when  H20  in 
the  plume  relaxes  to  its  ambient  values,  the  fate  of  liquid 
coverage  strongly  depends  on  ambient  conditions.  For 
example.  If  the  vapor  pressure  of  H20  on  alumina  particles 
exceeds  the  ambient  H2G  partial  pressure,  H20  may  evap¬ 
orate  changing  y  from  0,1-1  to  0,02. 

[19]  Perhaps,  the  most  surprising  results  of  the  Athena-2 
plume  measurements  is  a  persistence  of  condensed  H20  and 
HNO3  on  alumina  particles  under  temperature  conditions 
that  are  considerably  higher  than  the  NAT  or  ice  thresholds. 
While  the  evaporation  of  H20  from  alumina  particles  may 
be  decreased  by  the  presence  of  HNOf  [e.g,,  Warshawsky  et 
aL,  1999],  it  is  not  clear  whether  this  is  enough  to  explain 
the  measurements.  The  calculations  here  represent  constant 
values  of  temperature  and  pressure  (208. 1  K  at  70  hPa  and 
200.2  K  at  97  hPa)  and  are  valid  in  the  lower  stratosphere 
only.  Since  a  dependence  of  7,  <*i,  ot2,  ru  and  r2  on  number 
of  H20  and  HNO3  monolayers  on  alumina  particles  is 
unknown,  our  calculations  assume  constant  values  of  these 
parameters.  Further  laboratory  measurements  are  required  to 
infer  these  dependences,  incorporate  them  in  a  model,  and 
check  our  results. 

[20]  An  additional  caveat  is  that  the  properties  of  the 
emitted  alumina  particles  may  differ  from  those  used  in 
laboratory  studies.  For  example,  current  laboratory'  data  are 
available  for  hexagonal  a-Al203  [e.g,,  Molina  et  aL,  1997] 
and  it  is  not  obvious  that  they  are  valid  for  cubic  y-Al203, 
which  may  be  formed  in  solid-fueled  rocket  exhaust  and  is  a 
more  reactive  form  of  alumina  [Dai  et  aL,  1997]. 

[21]  Our  analysis  also  shows  that  small  particles  (i.e.,  d  < 
0.7  pm)  account  for  almost  2%  of  emitted  alumina  mass 
based  on  the  size  distribution  measured  by  the  NMASS, 
FCAS,  and  MASP  instruments.  It  is  very  important  to  know 
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the  fraction  of  alumina  mass  in  submicron  range  [Danilin  et 
al 2001b],  since  these  small  particles  have  a  long  residence 
time  in  the  stratosphere  and  are  available  longer  for  back¬ 
ground  chlorine  activation  via  (Rl)  on  their  surface,  thus 
affecting  ozone.  Assuming  that  Space  Shuttle  and  Titan-IV 
rockets  have  the  same  size  distribution  of  alumina  particles 
as  Athena-2  and  applying  results  from  the  global  model 
[Danilin  et  al,y  2001b,  Table  2],  the  ozone  depletion 
potential  (ODP)  of  alumina  emitted  by  these  rockets  could 
be  as  large  as  0.1  for  y  =  0.02.  Our  study  shows  that 
alumina  particles  may  be  covered  by  water  and  nitric  acid, 
raising  the  possibility  that  the  (Rl)  reaction  probability  y  is 
larger  than  0.02  under  background  conditions.  If  so,  the 
estimated  ODP  of  alumina  particles  may  exceed  the  value  of 
0.2.  Further  laboratory  and  field  measurements  are  required 
in  order  to  better  quantify  alumina  particle  effects  on  the 
ozone  layer. 

3,  Summary 

[22]  Athena-2  rocket  plume  measurements  during  the 
ACCENT  campaign  show  that  H20  and  HNO3  are  con¬ 
densed  onto  alumina  particles.  Quantifying  their  uptake 
using  the  AER  plume  model  constrained  by  all  relevant 
measurements,  we  found  that  (1)  the  H20  uptake  coefficient 
is  larger  than  3  x  10“4,  (2)  HN€)3  molecules  reside  for  5- 
52  min  on  A1203  surfaces,  and  (3)  alumina  particles  in  the 
background  atmosphere  may  be  covered  by  100-200 
monolayers  of  H20  and  0.1-10  monolayers  of  HN03. 
These  values  are  uncertain  by  at  least  a  factor  of  two  and 
require  additional  laboratory  studies.  Once  the  interactions 
of  alumina  particles  with  H20,  HNO3,  and  other  acids  are 
understood  based  on  laboratory  studies  and  field  measure¬ 
ments,  new  model  calculations  should  be  made  In  order  to 
estimate  effects  of  alumina  particles  on  the  atmosphere. 
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Observation  of  Stratospheric  Ozone  Depletion  Associated  With 
Delta  II  Rocket  Emissions 
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Abstract,  Ozone,  chlorine  monoxide,  methane,  and  submicron 
particulate  concentrations  were  measured  in  the  stratospheric 
plume  wake  of  a  Delta  II  rocket  powered  by  a  combination  of 
solid  (NH4CIO4/AI)  and  liquid  (LOX/kero^ene)  propulsion 
systems.  We  apply  a  simple  kinetics  model  describing  the 
main  features  of  gas-phase  chlorine  reactions  in  solid 
propellant  exhaust  plumes  to  derive  the  abundance  of  total 
reactive  chlorine  in  the  plume  and  estimate  the  associated 
cumulative  ozone  loss.  Measured  ozone  loss  during  two 
plume  encounters  (12  and  39  minutes  after  launch)  exceeded 
the  estimate  by  about  a  factor  of  about  two.  Insofar  as  only  the 
most  significant  gas-phase  chlorine  reactions  are  included  In 
the  calculation,  these  results  suggest  that  additional  plume 
wake  chemical  processes  or  emissions  other  than  reactive 
chlorine  from  the  Delta  II  propulsion  system  affect  ozone 
levels  in  the  plume. 

1.  Introduction 

Rockets  used  by  the  space  launch  industry  employ  a 
number  of  different  propellant  combinations  and  emit  a  variety 
of  exhaust  products  directly  into  the  stratosphere.  While 
rocket  exhaust  currently  represents  a  small  fraction  of  the  total 
impact  of  industrial  activity  on  the  stratosphere,  prudence 
requires  that  we  carefully  evaluate  its  impact.  Space  launch 
traffic  in  general  is  widely  expected  to  increase  during  coming 
decades  so  that  the  stratospheric  impact  of  rocket  emissions 
from  all  propellant  types  will  increase.  Previous 
investigations  have  focused  mainly  on  the  stratospheric 
impact  of  large  solid  rocket  motors  (SRMs)  using  ammonium 
perchlorate  {NH4CIO4)  oxidizer.  Recent  stratospheric 
measurements  have  shown  that  large  solid  fueled  rockets 
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dramatically  perturb  stratospheric  ozone  in  regions  covering 
up  to  several  hundred  square  kilometers  for  several  hours  after 
launch  [Ross  et  al.,  1997b], 

Solid  fueled  rockets,  however,  represent  only  a  portion  of 
the  emission  inventory  of  the  space  launch  industry.  Another 
widely  used  propellant  combination  is  liquid  oxygen  and 
kerosene  (LOX/RP;  RP  refers  to  RP-1  or  RG-1,  kerosene 
distillations  widely  used  as  rocket  fuel.)  Careful  study  of 
LOX/RP  combustion  emissions  Is  justified  for  several  reasons. 
LOX/RP  exhaust  accounts  for  a  significant  fraction  (about  one 
fourth  by  mass  in  1998)  of  the  total  stratospheric  emission  by 
rockets  and  several  powerful  LOX/RP  rockets  in  development 
will  increase  the  relative  impact  of  these  emissions.  In  this 
paper  we  present  measurements  obtained  in  the  stratospheric 
plume  wakes  of  a  Delta  II  rocket  that  provide  evidence  that 
reactive  chlorine  gas  may  not  be  the  only  chemically  active 
component  of  the  Delta  II  emission. 


2.  Plume  Wake  Measurements  and  Analysis 

The  Delta  II  plumes  were  sampled  following  daytime  launches 
on  November  7,  1996  (1600  UT)  from  Cape  Canaveral  Air 
Station  (28*  33’  N,  80*  18*  W)  and  on  May  17,  1998  (2316 
UT)  from  Vandenberg  Air  Force  Base  (120*  3T  W,  34*  48*  N). 
The  altitudes  of  the  plume  encounters  varied  between  18  and 
18.6  km.  Measurements  included  O3  [Sen  et  al,  1996] 
(University  of  Houston)  and  total  (volatile  and  nonvolatile) 
aerosol  concentrations  [Ross  et  al.,  1999]  in  1996  and  CIO 
[Pierson  et  al.,  1999],  CH4,  three  Independent  O3 
measurements  [Sen  et  al,  1996]  (University  of  Houston); 
[Proffitt  et  al.,  1989]  (NOAA);  and  [Proffitt  and  Rawlins, 

1998]  (PSI)  and  total  aerosol  concentrations  [Ross  et  al., 

1999]  in  1998.  Following  each  launch,  the  WB-57F  aircraft 
intercepted  the  plume  wake  six  times.  Only  the  1998  data  is 
comprehensive  enough  to  support  significant  conclusions 
regarding  the  details  of  plume  wake  chemistry  and  we  focus  on 
that  data  shown  in  Figure  1.  The  plume  encounters  are  clearly 
visible  as  large  departures  from  the  background  for  all  species 
except  €R$,  which  shows  a  modest  decrease  only  during  the 
fifth  encounter.  The  sixth  encounter  does  not  conform  with  the 
trends  in  composition  (progressively  greater  ozone  loss)  and 
duration  (progressively  longer)  observed  in  the  preceding 
encounters,  apparently  because  the  aircraft  Intercepted  an  edge 
of  the  main  plume.  Since  the  encounters  occurred  at  slightly 
different  altitudes  (between  18  and  18.6  km),  the  aircraft 
sampled  distinct  plume  segments. 

Representative  of  the  first  three  encounters.  Figure  2  shows 
expansions  of  the  second  encounter  when  the  plume  was  12 
minutes  old.  At  the  center  of  the  generally  symmetric  plume, 
ozone  concentration  was  reduced  by  70%  from  the 
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Figure  1.  Time  series  of  measurements  during  the  46  minute 
period  that  included  the  six  plume  encounters.  Units  are 
mixing  ratio  in  parts  per  million  by  volume  (ppmv),  parts  per 
billion  by  volume  (ppbv),  and  total  (volatile  and  nonvolatile) 
aerosol  number  density  in  the  size  range  0,3  to  4  pm.  The 
small  aerosol  event  after  the  fourth  plume  encounter  is 
interpreted  as  pi  encounter  with  the  WB-57F*s  own  exhaust. 
The  bottom  axis  is  Universal  Time  on  17  May  1998. 


background  value.  The  maximum  CIO  concentration  was  45  ± 
15  ppbv,  about  a  factor  20  greater  than  in  the  South  Polar 
springtime  vortex  [Pierson  et  a!.,  1999].  Methane 

concentration  in  the  plume  was  not  different  from  the  ambient 
background,  A  proposed  description  of  gas  phase  SRM  plume 
photochemistry  features  a  CIO  dimer  catalytic  cycle  unique  to 
warm  (  »  215  K)  SRM  plumes  [Ross,  1996;  Ross  et  al,  1997b]: 

2x  (Cl  +  O3  — »  CIO  +  O2  ) 

CIO  +  CIO  +  M  — >  CI2O2+  M 

Ci202+Cl-*Cl2  +  C100  (1) 

Cb  +  fty  2CI 
ClOO  +  M C1+ <V  M 
Net:  2O3  — ¥  30  2 . 

The  set  of  reactions  (1)  can  reasonably  be  taken  to  describe  the 
main  aspects  of  SRM  plume  kinetics  because  all  of  the  chlorine 
species  are  in  significant  excess  (  =  102)  of  ambient  species  that 
might  play  a  role  in  the  plume  kinetics  (HOx  and  NOx  for 
example)  and  because  all  other  Clx  family  reactions  are  slow 
compared  to  those  that  make  up  (1).  Our  desire  is  not  to 


explain  all  the  details  of  the  plume  kinetics  (nor  is  that 
presently  possible,  given  the  many  uncertainties  regarding  the 
actual  Delta  II  emissions).  Rather,  we  seek  to  evaluate  the 
ability  of  the  reaction  set  (1)  to  explain  the  observed  ozone 
loss. 

Within  a  few  minutes  after  launch,  the  principal  chlorine 
species  are  In  a  slowly  varying  steady-state  and  the  reactions 
(1)  can  be  solved  to  yield  a  relationship  between  [CIO]  (in 
units  of  cm'3)  and  total  reactive  chlorine  [Clx]  =  [CIO]  +  [Cl]  + 
2[C12]  +  2[C1202]: 

[Clx]  ~  1.5  [CIO]  +  5,4  x  ICT"  [CIO]2  (2) 

The  numerical  coefficients  were  determined  using  DeMore  et 
al.  [1997]  for  temperature  and  pressure  equal  to  215  K  and 
65  mbar,  respectively,  and  ambient  03  mixing  ratio  of  2  ppmv; 
the  coefficients  are  not  greatly  sensitive  to  these  conditions. 
Detailed  kinetics  calculations  [Rawlins  et  al.,  1998]  verify  the 
analytic  result  (2)  and  show  that  about  90%  of  the  Clx  is  in 
the  form  of  Cl2  during  the  steady-state.  We  average  the  plume 
ingress  and  egress  data  (to  reduce  the  influence  of  small  scale 
structure  in  the  plume)  and  fit  a  linear  expression  to  describe 
the  radial  dependence  of  [CIO]  and,  therefore,  {Clx]  using  (2). 
Integrating  derived  [Clx]  across  the  plume  and  taking  into 
account  the  known  uncertainty  in  measured  [CIO]  we  find  that 
the  total  abundance  of  reactive  chlorine  in  the  plume,  {Clx} 
(brackets  indicating  the  horizontal  Integration),  does  not 
exceed  0,9  moles  per  meter  of  altitude  (Mm4), 

When  Clx  is  greatly  In  excess  of  O3,  Cl2  photolysis  (J  = 

O, 22  min*1  for  the  relevant  local  conditions,  P,  F,  ZIttel  and 

P.  K.  Swaminathan,  personal  communications)  controls  the 
rate  of  ozone  destruction  and  the  horizontally  Integrated 
cumulative  ozone  loss  in  the  plume  and  {AO3}  during  a  time 
At  can  be  approximated  by 


Time  (10*  Seconds) 


Figure  2.  Expanded  time  series  for  the  second  plume 
encounter.  In  the  topmost  panel  data  from  the  NOAA  and  PSI 
instruments  are  given  by  circles  and  squares,  respectively. 
Data  rates  are  1  Hz  except  for  CEL  which  is  0.5  Hz.  The  bottom 
axis  Is  Universal  Time  on  17  May  1998,  WB-57F  airspeed  Is 
about  190  ms'1. 
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[i]  In  situ  measurements  of  NO,  NOy,  HN03,  CIO,  C02, 
H20  and  particles  were  made  in  the  plume  of  an  Athena  II 
solid-fiiel  rocket  motor  (SRM)  with  instruments  onboard  the 
NASA  WB-57F  high-altitude  aircraft.  By  normalization  to 
C02,  the  NOy  emission  index  was  calculated  to  be  2,7±0.6  g 
N02  (kg  fuel)-1,  HN03  is  a  significant  and  evolving 
component  plume  NOr  The  principal  source  of  HN03  is 
thought  to  be  the  heterogeneous  reaction  between  C10N02 
and  HC1  on  emitted  alumina  particles.  These  measurements 
provide  the  first  experimental  evidence  for  the  production  of 
HN03  and  the  activation  of  chlorine  on  SRM-emitted 
alumina  particles  in  the  lower  stratosphere.  INDEX 
TERMS:  0305  Atmospheric  Composition  and  Structure:  Aerosols 
and  particles  (0345,  4801);  0322  Constituent  sources  and  sinks; 
0340  Middle  atmosphere — composition  and  chemistry. 
Citation;  Popp,  P.  J.,  et  a!..  The  emission  and  chemistry  of 
reactive  nitrogen  species  in  the  plume  of  an  Athena  II  solid-fuel 
rocket  motor,  Geophys .  Res.  Lett 29(18),  1887,  doi:l  0.1029/ 
2002GL0 15197, 2002. 

1*  Introduction 

[2]  Current  SRMs  release  energy  through  the  oxidation 
of  aluminum  in  the  presence  of  ammonium  perchlorate 
(NH4CIO4),  resulting  in  the  emission  of  combustion  prod¬ 
ucts  directly  into  the  stratosphere  during  launch  [Prather  et 
ai ,  1990].  Nitric  oxide  (NO)  is  produced  directly  by  the 
SRM  combustion  process  and  in  secondary  combustion 
processes  known  as  afterburning  when  ambient  air  is 
entrained  into  the  hot  exhaust  plume  [Bennet  and  McDo¬ 
nald ,  1994],  NO  production  and  reaction  in  SRM  plumes  is 
not  well  studied.  Although  an  SRM  NO  emission  index  of 
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7  g  of  NO  (as  equivalent  N02)  per  kg  of  fuel  has  been 
reported,  this  value  did  not  account  for  the  post-emission 
oxidation  of  NO  in  the  plume  [Potter,  1 977], 

[3]  Observations  of  chlorine  monoxide  (CIO)  in  a  strato¬ 
spheric  SRM  plume  [Ross  et  ai,  2000]  suggest  that  emitted 
NO  is  oxidized  and  subsequently  forms  chlorine  nitrate 
(C10N02).  Laboratory  studies  have  shown  that  C10N02 
could  then  react  with  hydrogen  chloride  (HC1)  on  alumina 
(A1203)  particles  via 

AL03 

CIONO2  +  HC1  2L  Cl2  +  HNO3  (R1 ) 

[Molina  et  at. ,  1997]  to  effectively  convert  NO  to  nitric 
acid  (HN03),  It  should  also  be  noted  that  (Rl)  provides  a 
means  of  converting  HCI,  the  primary  chlorine  species 
emitted  by  SRMs  [Prather  et  ah,  1990],  to  active  chlorine. 
Hie  severe  ozone  (03)  loss  observed  locally  in  SRM 
plumes  is  known  to  result  from  reactions  involving  active 
chlorine  species  [J?e>5S  et  ai ,  1997a,  1997b].  While  it  has 
been  suggested  that  heterogeneous  chlorine  activation 
could  play  a  role  in  plume  03  loss  [ifcss  et  ai,  2000], 
the  contribution  of  (Rl)  to  active  chlorine  in  the  plume 
requires  further  observational  support. 

[4]  We  report  here  in  situ  measurements  of  total  reactive 
nitrogen  (NOy  =  NO  +  N02  +  C10N02  +  2N205  +  HN03 
. . .),  NO,  and  HN03  in  the  stratospheric  plume  of  an  Athena 
II  SRM.  No  such  in  situ  measurements  of  HN03  and  NOy 
have  been  reported  previously.  These  measurements  are 
used  to  calculate  an  NOy  emission  index  for  the  SRM  and 
to  infer  the  occurrence  of  (Rl)  on  SRM-emitted  alumina 
particles  in  the  lower  stratosphere. 

2,  Results  and  Discussion 

[5]  The  Athena  II  SRM  plume  was  sampled  following 
launch  from  Vandenberg  Air  Force  Base,  CA  (34°48'  N, 
120°37/  W)  on  24  September  1999  at  66060  s  Universal 
Time  (UT).  In  situ  measurements  were  conducted  onboard 
the  NASA  WB-57F  high-altitude  research  aircraft  as  part 
of  the  NASA  Atmospheric  Chemistry  of  Combustion 
Emissions  Near  the  Tropopause  (ACCENT)  mission.  Spe¬ 
cies  measured  include  NO  and  NOy  [Weinheimer  et  ah, 
1993],  HN03  [Neuman  et  al. ,  2000],  CIO  [/fas?  et  al , 
2000],  carbon  dioxide  (C02)  [Gates  et  al. ,  2002],  water 
vapor  (H20)  [May,  1998],  and  particles  [O.  Schmid  et  al, 
manuscript  in  preparation,  2002].  The  WB-57F  rocket 
plume  sampling  methodology  has  been  presented  in  detail 
elsewhere  [Rosy  et  al ,  1997a,  1997b].  The  plume  was 
sampled  a  total  of  6  times,  with  the  first  5  intercepts 
occurring  at  altitudes  between  18.3  and  18.5  km  in  the 
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lower  stratosphere.  The  final  plume  encounter  occurred  at 
16.5  km  during  the  WB-57F  final  descent,  and  is  not 
considered  in  any  further  analysis  here. 

2*1 .  NOy  Emission  Index 

[6]  Time  series  measurements  obtained  during  the  5 
plume  intercepts  are  shown  in  Figure  1.  Plume  ages,  as 
measured  from  the  time  of  emission,  ranged  from  3.8  to 
26.3  min.  The  duration  of  the  intercepts  varied  from  2.2  to 
12.9  s,  corresponding  to  flight  paths  of  0.4  to  2.1  km, 
respectively.  To  reach  these  dimensions,  an  SRM  exhaust 
plume  must  undergo  near-complete  dilution  with  ambient 
air.  NOy,  HN03,  C02,  and  CIO  mixing  ratios,  however,  still 
increased  substantially  above  background  values  in  the 
plume.  C02  values  are  shown  as  the  amount  above  the 
stratospheric  background  value  of  approximately  368  ppmv. 
NO  values  in  the  plume  as  soon  as  3.8  min  after  emission 
indicate  a  depletion  from  the  background  value  of  approx¬ 
imately  420  pptv  to  125  pptv  at  the  time  of  the  first  intercept 
and  to  smaller  amounts  in  the  remaining  Intercepts.  This 
measurement  illustrates  the  rapid  chemical  conversion  of 
NO  to  higher  oxides  of  nitrogen  in  the  plume  and  empha¬ 
sizes  the  importance  of  reporting  an  NOy  emission  index 
using  a  measurement  of  all  NOy  species. 

[7]  In  calculating  an  NOy  emission  index,  the  plume 
dilution  and  expansion  exhibited  in  Figure  1  is  accounted 
for  by  normalization  to  a  conserved  species  with  a  known 
emission  index  such  as  C02.  Since  the  SRM  emission 
index  of  C02  is  known  to  be  382  g  (kg  fuel)-1,  based  on 
the  carbon  content  of  the  fuel  burned,  measured  NOy  can 
be  directly  related  to  fuel  use  by  the  relative  abundances  of 
NOy  and  C02  in  the  plume.  Comparing  the  emitted 
quantities  of  NGy  and  C02  during  the  plume  encounters 
shown  in  Figure  1  is  accomplished  by  integrating  both 
species  after  subtracting  their  respective  background  val¬ 
ues.  The  integral  comparisons  are  meaningful  because 
each  in  situ  measurement  is  made  continuously  or  near- 
con  tinuously  in  a  constant  sample  flow  [Fahev  et  al, , 
1995a,  1995b]. 

[s]  The  plume-integrated  values  of  NOy  and  C02  (ANOy 
and  AC02,  respectively)  during  each  of  the  5  plume 
intercepts  are  shown  in  Table  1.  Although  NOy  and  C02 
are  conserved  in  the  plume  downstream  of  the  afterburning 
region,  the  Integrated  values  of  both  species  vary  by  more 
than  a  factor  of  2  during  the  5  intercepts.  This  variance 
results  from  both  the  differences  in  the  WB-57F  flight  paths 
within  the  plume,  and  the  non-uniform  dispersion  of  the 
plume  In  the  lower  stratosphere  [Belting,  2000].  Knowing 
ANOy  and  AC02  for  each  plume  intercept,  the  NOy 
emission  index,  EI(NOy),  can  be  calculated  according  to 

EI(N°y)  =  EI(C02)  •  *  ~  g  N03(kg  propellant)-1  (1) 

[Fahey  et  al ,  1995a]  where  EI(C02)  is  the  SRM  C02 
emission  index  and  (46/44)  represents  the  molar  mass  ratio 
of  NOy  (expressed  as  equivalent  mass  of  N02)  to  C02. 

[9]  Using  equation  1,  EI(NOy)  values  calculated  for  each 
of  the  5  plume  events  yield  an  average  of  2.7±0.6  g  N02  (kg 
fuel)-1  (Table  1).  This  value  is  higher  than  the  prediction  of 
1 .2  g  N02  (kg  fuel)” 1  calculated  by  one  of  the  authors  (RF.Z.) 
using  a  plume  flowfield  model  of  the  Athena  II  SRM 
initialized  for  conditions  at  18.6  km  altitude.  For  related 


Figure  1.  In  situ  measurements  of  NOy,  HNO3,  NO,  C02, 
and  CIO  in  the  plume  of  an  Athena  II  SRM  launched  from 
Vandenberg  Air  Force  Base,  CA  on  24  September  1999  at 
66060  UT.  The  panels  showing  data  during  the  first  plume 
intercept  at  66380  UT  represent  10  s  of  data.  Panels  for  the 
remaining  plume  intercepts  have  the  same  horizontal  scale. 
Plume  ages  at  the  time  of  Intercept  are  shown  at  the  top  of 
each  column. 

examples  of  SRM  plume  chemistry  calculations  using  super¬ 
sonic  SRM  exhaust  flowfield  models,  see  Zittel  [1994]  and 
Denison  et  al  [1994].  Approximately  30%  of  the  modeled 
NO  production  occurs  in  the  SRM  combustion  chamber 
where  the  conditions  affecting  NO  production  are  not  well 
understood.  The  remaining  NO  production  occurs  primarily 
via  the  Zeldovich  thermal  mechanism  [Blazowski  and  Saw¬ 
yer,  1975]  in  the  afterburning  region  of  the  plume  and  is 
highly  sensitive  to  local  temperature  and  mixing.  Because  of 
the  uncertainties  in  afterburning  plume  conditions  introduced 
by  approximations  to  supersonic  turbulent  mixing,  as  well  as 
the  possibility  that  additional  NO  production  paths  may  be 
discovered  to  be  significant  in  the  unusual  SRM  plume 
environment,  the  modeled  NO  production  Is  considered  an 
estimate  and  is  easily  subject  to  a  factor  of  2  uncertainty. 
While  the  model  does  not  follow  the  reactions  of  NO 
subsequent  to  afterburning,  the  calculated  NO  residue  can 
nonetheless  be  considered  the  sum  of  the  NOy  species  at  all 
later  times  In  the  plume  and  Is  therefore  directly  comparable 
to  the  experimentally  determined  value.  The  comparable 
EI(NOy)  values  from  the  measurements  and  model  predic- 
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Table  1,  Details  of  the  WB-57F  Intercepts  of  the  Athena  II  SRM  Plume 


Plume 

Int 

UT 

(s) 

Age3 

(min) 

T 

(K) 

P 

(hPa) 

TASb 

(mis) 

Atc 

(s) 

ANOyd 

(ppbv's) 

ACO,d 

(ppmvs) 

EI(NOy) 

(g  N02/kg) 

AHNG3d 

fppbv-s) 

HN03/NOy 

Observed 

HN03/NOy 
Lower  Limit 

1 

66380 

3.8 

208 

70 

202 

2.2 

178 

20.3 

3.5 

114 

0.64 

0.12 

2 

66733 

9.7 

208 

69 

191 

5.7 

125 

17.8 

2.8 

147 

1.18 

0.23 

3 

67077 

15.4 

207 

67 

173 

4.2 

80 

11.6 

2.8 

93 

1. 16 

0.22 

4 

67370 

20.3 

207 

69 

181 

8.7 

122 

24.6 

2.0 

152 

1.25 

0.24 

5 

67735 

26.3 

207 

66 

166 

12.9 

192 

31.6 

2.4 

243 

1.27 

0.24 

“Assumes  passage  of  Athena  II  SRM  through  intercept  altitudes  at  66154  UT. 
bTrue  air  speed  of  the  WB-57F  at  the  time  of  the  plume  intercept 

^Elapsed  time  of  the  plume  intercept,  defined  as  the  full-width  at  half-maximum  of  the  CIO  instrument  response  in  the  plume. 

A  indicates  plume-integrated  value  with  background  atmospheric  values  subtracted.  At  intercept  altitudes,  the  fuel  mass  flow  rate  of  the  Athena  II  SRM 
second  stage  was  590  kg  s  \  the  vertical  velocity  was  396  m  s"1  and  the  flight  path  angle  was  48°  from  vertical. 


tions  represent  satisfactory  agreement  considering  the  uncer¬ 
tainties  in  simulating  the  combustion  chamber  and  after¬ 
burning  plume  conditions  responsible  for  NO  production. 

2.2.  HN03  Production 

[10]  The  relative  plume-integrated  values  of  HN03 
(AHN03)  and  NOy  shown  in  Table  1  indicate  that  HN03 
Is  a  significant  and  evolving  component  of  NOy  in  the  aging 
rocket  plume.  The  amount  of  HN03  produced  within  the 
SRM  and  the  afterburning  plume  is  uncertain  and  is  not 
addressed  directly  by  the  plume  tlowfield  model.  Nonethe¬ 
less,  any  increase  in  HN03  observed  during  the  sampling 
passes  must  be  produced  in  situ  in  the  plume.  It  is  also 
possible  that  all  of  the  HN03  is  formed  in  this  manner.  Such 
production  could  occur  in  the  gas  phase  or  heterogeneously 
on  emitted  alumina  particles.  The  most  likely  gas-phase 
production  pathway, 

NO2+OH-+HNO3  (R2) 

is  expected  to  be  negligible,  since  it  requires  OH  to  be  in 
excess  of  I  ppbv  to  produce  nominal  HN03  values. 
Although  there  are  no  direct  measurements  of  OH  in  the 
Athena  II  plume,  the  flowfield  model  predicts  a  negligible 
residue  of  OH  in  the  cool  plume.  If  the  model  is  accurate, 
the  OH  density  is  too  low  by  several  orders  of  magnitude  to 
explain  the  observed  levels  of  HN03. 

[n]  The  heterogeneous  production  of  HN03  has  been 
demonstrated  In  laboratory  studies  to  occur  via  (Rl)  on  the 
surface  of  alumina  particles  [Molina  et  ah,  1997].  NO 
produced  directly  by  the  SRM  would  be  very  short  lived  In 
the  young  (<3.8  min  old)  plume  and  efficiently  converted  to 
C10N02.  Measurements  in  the  plume  indicate  that  less  than 
0.1%  of  the  emitted  NO  remains  after  3.8  min,  while  at  the 
same  time  NOy  levels  exceed  100  ppbv  (Figure  1 ).  This  rapid 
NO  removal  likely  occurs  by  reaction  with  CIO  in  the  plume: 

N0  +  CI0^N02+C1  (R3) 

With  15  ppbv  of  CIO  present  (the  approximate  mean 
abundance  of  CIO  in  the  plume  during  the  first  intercept,  see 
Figure  1),  the  lifetime  of  NO  in  (R3)  is  approximately  1  s. 
The  subsequent  formation  of  C10N02,  via 

M 

N02  +  CI0-C10N02  (R4) 

results  in  an  N02  lifetime  of  approximately  23  s.  It  should 
be  noted  that  the  lifetimes  of  NO  and  N02  stated  here  are 
upper  limit  estimates,  considering  that  the  CIO  abundance 
early  In  the  plume  was  likely  higher  than  the  measured  value 
of  1 5  ppbv  at  3.8  min.  The  net  result  of  (R3)  and  (R4)  is  that 
NOy  not  already  in  the  form  of  HN03  is  rapidly  sequestered 


into  the  relatively  stable  reservoir  species  C10N02  in  the 
young  SRM  plume.  Using  a  reaction  probability  of  0.02  for 
(Rl)  [Molina  et  al ,  1997],  and  an  observed  alumina  particle 
surface  area  density  in  the  plume  of  1.8 TO3  pm2  cm-3  (data 
not  shown),  the  lifetime  of  CI0N02  through  (Rl)  in  the 
young  plume  is  8.6  min.  Thus,  substantial  HN03  can  be 
produced  via  (Rl)  In  the  Athena  II  SRM  plume  by  the  time 
of  the  first  intercept  at  a  plume  age  of  3.8  min. 

[12]  The  observed  HN03  to  NOy  ratio  is  0.64  during  the 
first  plume  encounter  and  increases  above  unity  for  the 
remaining  intercepts  (Table  I).  Although  the  increases 
above  unity  are  within  the  combined  uncertainty  of  the 
measurements,  ratios  less  than  or  equal  to  one  are  expected 
based  on  measurements  made  prior  to  the  plume  intercepts 
and  during  a  previous  flight  [Neuman  et  al,  2000].  The 
HNO3  to  NOy  ratios  greater  than  unity  observed  here  are 
interpreted  to  be  the  result  of  some  fraction  of  the  HNO3 
being  adsorbed  on  alumina  particles  present  In  the  exhaust 
plume.  Size-resolved  measurements  of  particles  in  the 
plume,  assumed  to  be  composed  of  alumina,  indicate  a 
tri-modal  distribution  of  particle  surface  area  (SA)  in  the 
plume.  The  largest  particle  mode,  representing  65%  of  the 
surface  area  in  the  plume,  consists  of  particles  with  a  mean 
diameter  of  1.7  pm.  The  remaining  particle  surface  area  in 
the  plume  resides  in  smaller  modes  at  0.08  pm  and  0.16  pm 
[O.  Schmid  et  al,  manuscript  in  preparation,  2002]. 

[13]  The  HN03  instrument  samples  through  a  forward¬ 
facing,  sub-isokinetic  inlet,  which  causes  particles  larger 
than  approximately  0.1  pm  in  diameter  to  be  significantly 
enhanced  in  the  sample  flow  [e.g.  Northway  et  al,  2002].  In 
contrast,  little  to  no  particle  enhancement  is  expected  in  the 
side-facing  NOy  inlet.  Thus,  if  HN03  is  adsorbed  on 
alumina  particles,  the  observed  HN03  mixing  ratio  in  the 
plume  may  exceed  that  of  NOy  when  HNO3  constitutes  a 
significant  fraction  of  the  NOy  present.  A  lower  limit  for  the 
HN03  to  NOy  ratio  In  the  plume  is  derived  by  assuming  that 
all  HN03  in  the  plume  is  adsorbed  on  alumina  particles. 
Using  a  calculated  enhancement  function  for  the  HN03  inlet 
[e.g.  Northway  et  al ,  2002]  to  adjust  for  particle  enhance¬ 
ment,  the  lower  limit  ratio  ranges  from  0. 12  at  a  plume  age 
of  3.8  min  to  0.24  at  a  plume  age  of  26.3  min  (Table  1).  If 
any  HN03  Is  in  the  gas  phase  the  effect  of  particle  enhance¬ 
ment  will  be  reduced,  and  the  HN03  to  NOy  ratio  will  be 
greater  than  the  lower  limit  values  and  approach  the 
observed  ratios.  The  actual  HN03  to  NOy  ratio  can  only 
approach  unity  if  almost  all  HN03  in  the  plume  Is  in  the  gas 
phase  and  the  effect  of  particle  enhancement  is  minimized. 

[14]  The  increase  in  the  HNO3  to  NOy  ratio  observed  in 
the  plume  between  the  first  and  second  plume  encounters  is 
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likely  due  to  continued  In  situ  production  of  HN03. 
Although  there  is  uncertainty  in  the  effective  surface  area 
density  in  the  plume,  the  lifetime  of  C10N02  In  the 
production  of  HN03  is  long  enough  that  (Rl)  could  still 
be  occurring.  The  near-constant  HN03  to  NOy  ratios  after 
the  second  encounter  are  likely  due  to  the  cessation  of 
HN03  production  in  the  plume.  Cessation  could  occur  If 
CIONO2  is  completely  converted  to  HN03  or  if  (Rl)  slows 
due  to  contamination  of  the  alumina  surface  by  adsorbed 
HN03,  as  suggested  by  Molina  et  al.  [1997]. 

[15]  As  a  final  consideration,  we  note  that  nitric  acid 
trihydrate  (NAT;  HNO30H2O)  is  not  stable  in  the  aging 
plume.  Water  vapor  mixing  ratios  are  near  20  ppmv  during 
the  first  plume  intercept  and  11-15  ppmv  for  the  later  cross¬ 
ings  (data  not  shown).  For  20  ppmv  H2O  and  the  measured 
plume  temperature  of  208  K  at  70  hPa,  HNO3  mixing  ratios 
must  exceed  approximately  140  ppbv  for  NAT  to  be  stable 
[Hanson  and  Mauersberger ,  1988].  This  value  far  exceeds 
the  maximum  HN03  mixing  ratio  of  5 1  ppbv  observed  during 
the  first  plume  intercept.  As  H20  mixing  ratios  drop  to  11  —  15 
ppbv  in  the  later  plume  crossings,  NAT  continues  to  be 
unstable  at  the  observed  HNO3  mixing  ratios.  Before  the  first 
crossing  (at  a  plume  age  of  3.8  min),  NAT  may  have  been 
thermodynamically  stable  for  a  brief  period  in  the  expanding 
plume  because  the  plume  cools  to  ambient  temperatures  In  a 
few  seconds  and  the  HN03  saturation  mixing  ratio  depends 
non-linearly  on  the  H20  mixing  ratio.  In  addition,  NAT  would 
be  present  at  all  the  plume  crossings  if  the  ambient  temper¬ 
atures  were  slightly  lower.  Comprehensive  plume  models 
should  include  the  possibility  of  NAT  formation  and  the 
associated  heterogeneous  reactions  on  NAT  particles. 

3.  Implications 

[i6j  The  results  presented  here  provide  the  first  measure¬ 
ments  of  the  NOy  emission  index  for  an  SRM,  highlighting 
the  importance  of  validating  the  parameterizations  and 
assumptions  used  in  afterburning  models.  The  accumulation 
of  NOy  emitted  by  SRMs  could  potentially  affect  03  pro¬ 
duction  and  loss  processes  in  the  stratosphere.  A  typical  SRM 
launch  scenario  described  by  Danilin  et  al,  [2001]  results  in 
the  combustion  of  3.4- 106  kg  of  fuel  in  the  stratosphere 
annually.  The  calculated  NOy  emission  index  of  2.7±0.6  g 
NO2  (kg  fuel)-1  presented  here,  therefore,  results  in  the 
annual  emission  of  9.2* ID3  kg  of  NOy  directly  into  the 
stratosphere.  In  contrast,  the  production  of  reactive  nitrogen 
in  the  stratosphere  from  the  decomposition  of  N20  is  approx¬ 
imately  1.5-1G9  kg  N02  yr-1  [Murphy  and  Fahey ;  1994]. 
Thus,  the  contribution  of  SRM  emissions  to  the  stratospheric 
NOy  budget  Is  negligible  compared  to  that  from  natural 
sources.  Future  NOy  emissions  from  other  rocket  propellant 
types  (i.e.  hydrazine-based  fuels)  may  not  be  negligible. 

[1  ?]  These  data  also  provide  the  first  experimental  evidence 
forthe  production  ofHN03  via  (Rl) on  SRM-emitted alumina 
particles  In  the  lower  stratosphere.  In  addition,  they  provide 
evidence  that  HN03  may  remain  adsorbed  on  the  alumina 
particles.  Although  the  occurrence  of  (Rl)  results  in  chlorine 
activation  in  the  plume,  the  importance  for  local  03  loss  is  not 
clear  since  active  chlorine  Is  otherwise  abundant  in  the  plume 
[Rays  et  al ,  1997a,  1997b].  With  the  potential  future  increase 
in  launch  vehicle  traffic,  understanding  how  reactions  on 
alumina  particles  affect  global  stratospheric  03  will  be  critical 
in  assessing  the  effects  of  the  space  launch  industry. 
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